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INFLUENCE OF VIBRATION MODES 
ON CONTROL SYSTEM STABILIZATION 
FOR 

SPACE -SHUTTLE TYPE VEHICLES 


By H. G. Greiner 


1 SUMMARY 


An investigation was made to determine the feasibility of using conventional auto- 
pilot techniques to stabilize the vibration modes at the liftoff flight condition for two 
different space shuttle configurations. The objective was to determine if these con- 
ventional techniques would still be adequate to stabilize the elastic modes for these 
configurations, which have greater modal density and asymmetry than any previous 
launch vehicle. One configuration is called the dual flyback vehicle in which both the 
orbiter and booster vehicles have wings and complete flyback capability. The other 
configuration is called the solid motor vehicle in which only the orbiter has flyback 
capability. This latter vehicle consists of four elements: two strap-on solid motors, 
an external liquid propellant tank and the winged orbiter. The results of the linear 
stability analyses for each of these vehicles are summarized below. 


Dual Flyback Vehicle 

Adequate stability margins were provided for all parasitic modes (propellant slosh- 
ing and elastic vehicle) at the liftoff flight condition. Considerable latitude was available 
to establish satisfactory control mode (rigid body) response characteristics. The propel- 
lant sloshing modes exhibited small positive damping, but indications were that they would 
become less stable with increasing flight time. Essentially their stability properties were 
similar to those that have been witnessed on current operational launch vehicles. 

Satisfactory elastic- mode stability margins were obtained for the symmetric mode 
case when the autopilot pitch control plane was configured with the rate sensor located 
aft of the first body mode antinodes. The forward rate sensor configuration, including 
the autopilot compensation required to stabilize the first body modes , exhibited one 
slightly unstable elastic mode. However, indications were that with careful attention 
to actuator system design, the required stability margins could be provided for this 
mode . 

A two -rate -sensor (forward and aft) configuration also resulted in satisfactory 
stability margins for all of the elastic modes. This concept would permit a dominant 



aft rate sensor configuration early in flight, switching to a dominant forward rate 
sensor configuration for later flight times when the modal stability with the aft rate 
sensor typically becomes marginal. 

In the antisymmetric mode case , the autopilot yaw control plane exhibited satis™ 
factory stability margins for both the aft rate sensor configuration and the forward 
rate sensor configuration, with the appropriate first body mode stability compensation 
included in the latter. The roll control plane was also satisfactorily configured with 
either a forward or aft mounted rate sensor. Coupling between the yaw and roll auto- 
pilot channels was small , with the result that elastic modes in the coupled system 
exhibited adequate stability margins. 

The autopilot configurations that provided satisfactory stability margins were not 
found to be unusually sensitive to variations in the elastic mode natural frequencies . 
Only small adjustments in the stability compensation were required to cope with these 
frequency changes. 

It was found that local structural deformations at the autopilot sensor locations , 
particularly the rate sensor, could significantly alter the stability characteristics that 
were obtained from equivalent centerline type modal data. This condition was briefly 
examined in the roll autopilot control plane using the three-dimensional booster anti- 
symmetric mode shapes. 

Two anticipated stability problems at the liftoff flight time did not materialize for 
this vehicle . Negligible coupling was evidenced between the control mode , the propel- 
lant sloshing modes, and the lowest frequency elastic modes. Also, the elastic mode 
frequency density did not appear to be as important a factor as originally anticiplated. 


Solid Motor Vehicle 

With the exception of the yaw/roll coupled system, positive stability was demon- 
strated for all parasitic modes at the liftoff flight condition. However, the stability 
margins for some of the elastic modes were not satisfactory. 


The relatively small control thrust available from the orbiter engines required 
high autopilot gains to obtain satisfactory control mode (rigid body) response character- 
istics . Propellant sloshing stability was not a factor in this analysis because of the 
very small sloshing masses relative to the total vehicle mass (less than 1 percent for 
the oxidizer and 1/2 percent for the fuel). The most dominant sloshing mode was the 
oxidizer mode , and it assumed a stable root locus configuration that was hot 
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significantly influenced by any of the other oscillatory modes or by the autopilot elastic 
mode stability compensation that was investigated. 

The mode shapes exhibited a high degree of interaction between the vehicle ele- 
ments through the links that join these elements together. The modal data from these 

mode shapes, together with the high autopilot gains, produced stability compensation 
requirements which, in many cases, were not practical with conventional autopilot 
techniques. 

For the symmetric mode case the first orbiter body modes (elastic modes exhibit- 
ing classic first mode shape on orbiter) could not be stabilized for either forward- or 
aft-mounted rate sensors on the orbiter vehicle. The extremely high modal gain in the 
forward rate sensor case was incompatible with conventional compensation techniques. 
The aft rate sensor case showed some promise since it exhibited a negative damping 
ratio of only 2 percent. However, the high modal gain and the opposing nature of the 
loci from two adjacent open-loop poles made stability compensation relatively ineffec- 
tive. A brief attempt to stabilize this aft rate sensor configuration by selected varia- 
tion of mode frequencies produced only a marginally stable condition. Finally the two- 
rate-sensor configuration was employed, and the results indicated that all elastic 
modes could be stabilized. This configuration exhibited the greatest potential for pro- 
viding adequate stability margins for these elastic modes. The stabilizing of these 
first body modes was significantly complicated by the location of the "tail-wags-dog" 
frequency in the midst of these first-mode frequencies. 

Autopilot sensor locations on the external tank and the solid rocket motors were 
also investigated. The results exhibited unstable modes that were not amenable to 
conventional compensation techniques in all cases. 

The effect of increasing the control forces was briefly investigated by adding a 
thrust vector control system to the solid motor rockets. The resulting reduction in 
the autopilot forward loop gain significantly reduced the elastic modal gains and 
improved the system stability for the case where the autopilot sensors were located 
on the orbiter. However, the stability complication posed by the location of the "tail- 
wags-dog" zero was not appreciably reduced. 

In the antisymmetric mode case , all elastic modes were stable for the autopilot 
yaw control plane with an aft rate sensor located on either the orbiter or the external 
tank. For the autopilot roll control plane, stable elastic modes resulted with an aft 
attitude sensor located on either the orbiter or the external tank and an aft rate sensor 
located on the external tank. However, for considerations similar to those encountered 
in the symmetric mode case, none of these stable configurations provided 
satisfactory margins of stability. 
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Although the stable planar autopilot configurations were used, the yaw /roll 
coupled system exhibited unstable elastic mode properties. The analysis results in- 
dicated that there was significant coupling through the yaw attitude control system and 
that this coupling was detrimental to system stability. 
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2 INTRODUCTION 


The objective of this study was to conduct a preliminary investigation of the 
feasibility of using conventional control system techniques to stabilize the vibration 
modes at the liftoff flight condition for space -shuttle type vehicles. Two quite different 
space shuttle concepts were considered, both having greater modal density and asym- 
metry than any previous launch vehicle. They were: 1) a piggy-back configuration with 
both the booster and the orbiter vehicles having complete fly -back capability, and 
2) a strap-on solid motor configuration in which only the orbiter vehicle has fly -back 
capability. In this report the first configuration is referred to as the dual flyback 
vehicle , and the second is called the solid motor vehicle . 

The study was divided into three parts. Part one derived the perturb ational 
equations that were used in the stability analyses. Part two generated a consistent 
set of basic data for each of the study configurations . Part three performed the 
"time slice" stability analyses for each of the vehicles at the liftoff flight condition. 

The stability analyses for each vehicle were divided into two parts , which reflected 
the constraints that were applied to the elastic modal data. The first part considered 
elastic modal data subject to symmetric vehicle boundary conditions, and the second 
part used elastic modal data subject to antisymmetric vehicle boundary conditions. 

In terms of the vehicle attitude control system, the first part included the pitch control 
plane, and the second part included the yaw and roll control planes. Where practical, 
variations in parasitic mode (propellant sloshing and elastic vehicle) natural frequencies 
and damping ratios were considered for their effect on the conventional autopilot 
parameters . 

The basic data used in the dual flyback vehicle analyses were obtained from the 
final results of an earlier space shuttle study. The solid motor vehicle basic data 
were primarily obtained from a follow-on to this study. However, due to the timing 
between the progress of this autopilot study and the advent of the solid motor vehicle 
concept, these data are not as refined as the data applicable to the dual flyback vehi- 
cle. This was particularly true for the critical elastic vehicle models. Thus, it was 
not considered appropriate at this time to conduct as detailed an analysis for this 
vehicle as was conducted for the dual flyback vehicle. Instead, the solid motor vehi- 
cle analysis attempted to illuminate the most important considerations in applying 
conventional control system techniques to this vehicle. 

As stated earlier, this study considered two vehicle configurations representing 
two quite different space shuttle concepts. The first configuration was the piggy-back 
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delta wing vehicle shown in Figure 1. Both the booster and orbiter vehicles have 
independent flyback capability. During the ascent flight phase , all forward thrust 
and all control is accomplished with the booster. Forward thrust is obtained from the 
booster's twelve 2,446, 000N (550, 000 lb) thrust liquid oxygen/liquid hydrogen rocket 
engines . 



WEIGHT 


(lb) 

BOOSTER-LIFTOFF 

1,900,000 

(4,188,000) 

ASCENT PROPELLANT 

1,534,000 

(3,382,000) 

CRUISE PROPELLANT (JP) 

65,300 

(144,000) 

LANDING 

298,500 

(658,000) 

AREA 

m 2 

(It 2 ) 

EXPOSED WING 

469 

(5,047) 

EXPOSED CANARD 

47 

(504) 

VERTICAL TAIL 

139 

(1,500) 



Figure 1 . Dual Flyback Vehicle 


The size of the booster body is primarily a function of the amount of liquid propel- 
lants required to satisfy the mission profiles. This particular configuration contains 
two 10 M (33 ft) diameter tanks. The forward tank contains the liquid oxygen, and 
it has a capacity of 1160 M 3 (41 , 000 ft 3 ). The liquid hydrogen is contained in the aft 
tank, which has a volume of 3400 M 3 (120,200 ft 3 ). 


During the ascent flight phase vehicle control is accomplished by gimballing the 
booster rocket engines. All twelve engines are gimballed in both the pitch and yaw 
control planes. For roll control all twelve engines are differentially gimballed in 
both pitch and yaw. This maximum roll control thrust capability is required by the 
high degree of aerodynamic yaw/roll coupling exhibited by these piggy -back con- 
figurations . 
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The second configuration was the piggy-back solid motor vehicle shown in Figure 
2. For this vehicle, only the orbiter lias independent flyback capability. Forward 
thrust is obtained from the two 3.96 M (156 in.) solid rocket motors, each providing 
12,450,000 N (2,800,000 lb) of thrust and the three 1,650,000 N (370,000 lb) thrust 
liquid propellant main orbiter engines. 



Figure 2 . Solid Motor Vehicle 

The orbiter external propellant tank contains two 7.62 M (25 ft) diameter tanks. 
The forward liquid oxygen tank has a volume of 462 M 3 (16,300 ft 3 ), and the aft liquid 
hydrogen tank volume is 1320 M 3 (46,500 ft 3 ). 

The study mainly concentrated on using only the three main orbiter engines for 
vehicle control. This is certainly the desired approach from both cost and weight 
considerations. All three engines are gimballed for pitch and yaw control, but only 
the two engines that are symmetrically located relative to the vehicle's pitch plane 
are differentially gimballed in pitch for roll control. 

Another control configuration for this vehicle uses thrust vector control of the 
solid motors in conjunction with the main orbiter engines. This study considered 
secondary injection TVC for the solid motors. Aerodynamic surfaces were not con- 
sidered because of their ineffectiveness at the liftoff flight condition. 
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3 ANALYSIS PROCEDURE 


This study constitutes a linear analysis of the influence of the launch vehicle vibra- 
tion modes on overall vehicle stability. Nonlinearities and nonlinear effects are not in- 
cluded in the dynamic models, nor in the interpretations of the study results unless 
specifically stated to the contrary. 

The first part of this section describes a conventional autopilot and the particular 
model that is used in this study. The major features of the dynamic models are dis- 
cussed briefly in the second part. The third part contains a discussion of conventional 
autopilot design considerations with particular emphasis on parasitic mode stability. 
This part presents the background for the results discussed in the next two sections. 
The discussion in the third part and the results of the study are presented in terms of 
the root locus analysis method. A brief description of this technique is contained in 
Appendix E . 


Conventional Autopilot 


The word "conventional" refers to non-adaptive autopilots having proportional 
plus derivative feedback. These control systems are comprised of independent con- 
trol channels, one each for the pitch, yaw, and roll control axes. A functional con- 
trol channel for a conventional autopilot is shown schematically in Figure 3 . 

The primary feedback loop for this vehicle attitude control system is through the 
attitude sensor (displacement gyro or platform). The attitude sensor produces a sig- 
nal proportional to the angular difference between the instantaneous vehicle attitude 
and the desired vehicle attitude. This error signal is then used to displace the thrust 
chambers and cause the vehicle to change attitude in a direction that nulls the error 
signal. Vehicle guidance steering is effected by changing the desired vehicle attitude. 

To obtain the required rigid body stability and response characteristics, lead 
compensation is incorporated through a rate feedback loop. The signals proportional 
to the vehicle rotational rates are provided by the rate sensors, one for each control 
axis. 


9 



ATTITUDE GAIN 


FORWARD LOOP STATIC GAIN 


k SF 


VEHICLE 

ANGULAR 

DISPLACEMENT 

DESIRED 
VEHICLE 
ATTITUDE 


ATTITUDE ! 
+Q SENSOR 


r 


VEHICLE 

ANGULAR 

RATE 




_ . _ _ 

I 

FEEDBACK 

TRANSDUCER 

1 

1! 


SIGNAL 

! J 

L- 

SERVO 

CONDITIONING 

hfO 

ACTUATION 

FUNCTION 

1 + 
1 


SYSTEM 


L. 


ENGINE 

POSITIONING 

SYSTEM 


ENGINE 

-►ANGULAR 

DISPLACEMENT 


. J 


RATE GAIN 


Figure 3 . Conventional Autopilot Schematic 

The total system error signal is the sum of the attitude reference error signal and 
the rate signal. This total error signal is sent through a conditioning function before 
being applied to the engine positioning system. The primary purpose of the condition- 
ing function is to provide additional stability compensation, as may be required, to 
prevent the buildup of oscillations in the autopilot loop that would otherwise occur from 
the coupling of parasitic modes into the control system rate and attitude sensors. 

Another purpose is to provide a signal that is approximately proportional to the time 
integral of the attitude reference error signal. This integral signal maintains the atti- 
tude error signal at null in the presence of steady-state torques (i.e. , eliminates the 
need for the vehicle to be "tilted" to compensate these steady-state torques). However, 
this function can be satisfactorily supplied by guidance in some instances. Part or all of 
this conditioning function could also be accomplished in each of the feedback loops for- 
ward of the sensor summing junction. 

The conditioned error signal is then fed to the engine positioning system, a second- 
ary position feedback control system. This system causes the actual angular displace- 
ment of the vehicle's control thrust chambers. 
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Figure 4. Attitude Control System Block Diagram 

Forward loop gain (end-to-end attitude gain): 

K - k k ra -d» engine displacement 
A D SF rad. attitude sensor input 

End-to-end rate gain: 


K A K R ^ k SF rad 


rad, engine displacement 
— rate sensor input 


K can now be determined in terms of the functional schematic gains as: 
R 

K A K R k SF k R 


R K, 


*D*SF k I 
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Thus, Kj^ is tiie rate “to -attitude gain ratio. This technique for expressing the block 
diagram gains enables the system rate gain dependence to be handled as an equivalent 
lead filter. 

The conditioning function contains stabilization filters whose principal function is 
to stabilize the higher frequency parasitic modes; viz. , the elastic (bending) modes. 

A secondary function is to attenuate high frequency signals; this function aids in noise 
suppression. The integration function is not included for the purposes of this study, 
since it will have little effect on the parasitic modes. 

The rate sensor (rate gyro) dynamic properties can be modeled by a quadratic lag 
function. The natural frequency for a typical rate gyro is 160 rad/sec with a damping- 
ratio of greater than 0.4. These dynamic properties will have an effect on only the 
higher frequency elastic modes , and they have not been included in the autopilot model. 
However, if they prove to be significant, these effects can be considered in evaluating 
the analysis results. The attitude sensor dynamic properties are not significant in the 
range of frequencies considered in this study. 


Analysis Models 


This section discusses the general features of the mathematical models used in 
this analysis. A detailed description of each of these models and the complete equa- 
tions of motion are contained in the appendices . 

For stability studies, the dependent variables are small perturbational quantities. 
The measured variables can be considered to consist of two parts, a steady-state com- 
ponent that varies slowly as a function of time and a perturbational component that 
varies relatively rapidly at frequencies corresponding to the various modes of oscilla- 
tion of the launch vehicle. These latter quantities are the subject of interest in dynam- 
ic stability studies. 

Since the symmetry of these launch vehicles corresponds to that of an airplane, 
the equations of motion follow closely the standard development used for aircraft 
dynamic analysis. The longitudinal set of equations describes the motion in the plane 
of symmetry, and the lateral equations describe the motion out of this plane. In the 
stability analysis these two sets of equations are considered to be uncoupled from each 
other. The equations are linearized by neglecting the products of perturbation variables, 
using the small-angle approximation for trigonometric functions, and setting the Euler 
angle rates equal to the body rates. 
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Elastic vehicle model. - Since the vehicle elastic modes are orthogonal, coup- 
ling appears only through the external forces. In the elastic models the engine masses 
have been lumped into the mass of the thrust structure. Thus, the only engine inertia 
terms included in the equations of motion are those associated with the engine gimbal 
angle variable. For all of the elastic modes, the damping ratio used in the equations 
is a relatively conservative value of one-half per cent. 

The distributed aerodynamic properties are not included in the bending mode 
generalized forces. The small value of dynamic pressure at liftoff indicates that 
these aerodynamic forces are negligible relative to the TVC system and engine 
inertia forces. Generally, the effect of the aerodynamic terms on the bending mode 
stability is to slightly increase the effective modal damping and decrease the mode 
frequency. However, they may be destabilizing to "interconnect" modes in a piggy- 
back configuration at later flight times. 

Two sets of modes were generated for each vehicle, reflecting the vehicle 
symmetry. The symmetric set portrays the elastic vehicle characteristics in the 
plane of symmetry, and the antisymmetric set describes these characteristics out 
of the plane of symmetry. 

Propellant sloshing model. - The spring-mass and pendulum analogies are 
widely used to describe the dynamic characteristics of propellant sloshing. For 
present operational launch vehicles the "mechanical" parameters that are used in 
these models are derived from the solution of the hydrodynamic equations. This 
solution assumes that the tank is axisymmetric in shape, the vehicle center of mass 
lies on the tank centerline , and the vehicle acceleration vector is aligned with this 
centerline. It is these latter two assumptions which are not directly applicable to 
the space shuttle configurations. 

A recent study (Ref. 1) investigated the effects of acceleration vector "tilt" re- 
lative to the tank centerline under one-g conditions. The results of this study indicate 
that the axisymmetric model parameters are valid in cylindrical tanks for "tilt" 
angles of approximately 0 . 087 rad (5 deg) and fluid levels with a greater than one 
height- to-diameter ratio. These conditions approximate the liftoff conditions for 
the shuttle vehicles . 

Based upon these results and the relatively small center- of-mass offsets at lift- 
off, the axisymmetric propellant slosh models are considered to be adequate at this 
time of flight. The pendulum analogy is used in these stability analyses for the dual fly- 
back vehicle booster propellant tanks and for the solid motor vehicle orbiter external 
tanks. Propellant sloshing is not considered in the dual flyback vehicle orbiter 
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propellant tallies because these tanks are "full" throughout the ascent flight phase 
(i.e. , the slosh masses are negligible). 

Only the first lateral slosh modes are considered in this study. This mode has 
the lowest frequency, the least inherent fluid damping, and the greatest effective mass 
of moving fluid. The masses associated with the second and higher slosh modes for 
these shuttle configuration tanks are only a few per cent of the first mode mass for all 
fluid levels except the near empty condition. 

Thrust vector control system models . — Equivalent lumped parameter engines and 
actuators are used to generate the pitch and yaw control forces. These equivalent en- 
gines are located on the appropriate vehicle centerline. Roll control forces for the 
dual flyback vehicle are modeled by an equivalent force couple centered about the 
booster centerline. This force couple is generated by multiplying the total booster 
thrust by a derived moment arm such that the control moment is equivalent to that ob- 
tained from gimballing all twelve booster engines differentially in both pitch and yaw. 
Since roll control forces for the solid motor vehicle are obtained by differentially 
gimballing two orbiter engines in pitch, this model is straightforward. 

A third-order actuator model is used in all three control channels. This model 
has a closed loop resonant frequency, at load, of 50 rad/sec and a damping ratio of 
three tenths. The servo resonant frequency is 70 radians per second and the servo 
loop gain is 13 sec -1 . These model parameters are based on the requirements that 
were generated in earlier space shuttle studies. 


Conventional Autopilot Design Considerations 


The control function of the flight control system provides acceptable stability and 
control throughout flight of the launch vehicle. In this respect, the terms stability and 
control take on different connotations. Control implies the ability of the vehicle to 
execute steering commands generated by guidance and programmed steering and, 
having executed these commands, to maintain the desired trajectory in the presence of 
external disturbances. Stability, on the other hand, relates to the stabilization of the 
"parasitic" modes. The "parasitic" modes, e.g. , the propellant sloshing and the 
elastic modes, arise from the many resonances caused by the physical properties of 
the launch vehicle. The stability function of the control system is to prevent parasitic 
modal oscillations from passing through the control channel with amplification and 
phasing that would cause the oscillation to be sustained or to diverge. It is this stabil- 
ity function that is of primary concern in this study. 
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The autopilot gains, and Kr, are primarily determined by the desired control 
mode parameters: rigid body frequency and damping ratio. The desired rigid body 
control response is dominated by the need for minimal (within practicable limits) 
dynamic overshoot in response to atmospheric disturbances. It is axiomatic in control 
theory that minimal response to disturbances is achieved with high gains and bandwidth 
in the feedback and compensation signal paths. Unfortunately, this is incompatible with 
the lightly damped structural modes and an acceptable compromise must be worked out. 

The control mode response is considered in this study only for the purpose of ob- 
taining "reasonable " autopilot gains with which to evaluate the stability of the para- 
sitic modes. No effort is made to optimize this response. Since increasing gain, 
specifically the product of and K^, is generally destabilizing to these parasitic 
modes, the smallest acceptable values are of particular interest. For this purpose a 
rigid body response with a damped frequency of approximately 1 to 1.5 rad/sec and a 
damping ratio of between 0. 3 and 0.4 is used as a guide. It should be kept in mind, 
though, that it may be necessary to increase these gains to obtain a "tighter" response 
during i the higher dynamic pressure region of flight. 

The principal concern during the liftoff flight condition is the stability of the 
higher-frequency parasitic modes associated with airframe flexibility. These modes 
of oscillation manifest themselves as a series of lightly damped poles and zeros in the 
transfer function representing vehicle dynamics. 

The response of the launch vehicle at high frequencies is usually dominated by the 
lowest-frequency flexible modes as these have the least inherent damping, the greatest 
gain through the flight control system, and are the closest in frequency to the atmos- 
pheric disturbances that are the potential source of their excitation. Vehicle vibration 
in one or more of these modes contributes to the bending moment experienced by the 
vehicle and constitutes a significant contribution to the total vehicle loads, especially 
for the largest class of launch vehicles where the excitation is greatest. For these 
reasons the stabilization of these modes is usually achieved by ensuring that the phase 
of the parasitic signal picked up by the system sensors will not cause regenerative 
oscillations. In other words, the locus emanating from the elastic mode pole is forced 
to depart into the left-half plane. It is even possible to phase this signal such that the 
thrust vector deflection acts to stabilize the oscillation; the closed -loop flexible mode 
root will have greater damping than the open-loop pole, having a beneficial effect on the 
vehicle, loads. 

The higher-frequency flexible modes, by contrast, have greater inherent damping, 
contribute less to the vehicle loads, and have less gain through the flight control system 
loop. These modes are gain stabilized as opposed to the phase stabilization of the lower 
modes; the bending oscillations are decoupled from the flight control system by high 
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attenuation at these higher bending frequencies. For these modes the closed-loop root 
will always be in the left-half plane for all phase orientations. In addition, the pre- 
cision of the modal data at the higher modal frequencies makes reliance on phase stabil- 
ization speculative; the uncertainties in the mathematical model are too great. 

Those flexible modes in the range of the engine position servo resonance are a 
special case. There may be sufficient friction in the engine gimbal to ensure satis- 
factory stability. If not, then these modes must be examined closely with the possibility 
in mind of modifying the position servo transfer function itself to obtain the required 
degree of stability. In particular the load torques become quite important. 

An important elastic mode stability consideration is the frequency of the so-called 
"tail -wags -dog" zero (TWD). This is the frequency at which the lateral force at the 
gimbal axis due to thrust forces is exactly balanced by the inertial force of thrust 
chamber motion. Below this frequency the thrust forces are dominant and above it the 
engine inertia forces are dominant.* The significance of this zero for elastic stability 
is that the locus departure angle from an elastic pole undergoes a 180-degree phase 
shift when the pole passes through this frequency. This makes phase stabilization of a 
mode whose open-loop pole frequency is very near the TWD zero highly impractical. 

The gain and phase of an elastic mode locus is strongly dependent upon the modal 
slopes at the attitude and rate sensor locations. Therefore, the choice of sensor loca- 
tion will have considerable influence on the stability properties of the flight control sys- 
tem. This situation will be illustrated by considering only the first body bending mode. 
The pole -zero plot for this configuration is shown in Figure 5. The first diagram 
shows the effect of rate and attitude sensors (located at the same point) being slightly 
aft of the first mode antinode. For this position the modal slope at the sensor location 
will have the same sign as the slope at the engine gimbal point. In the second diagram 
the sensors are slightly forward of the antinode and the sensor modal slope has the 
opposite sign to the slope at the engine gimbal point. Root locus departure angles (the 
direction in which the root locus emanates from the bending pole) are shown assuming 
flat (zero phase shift) autopilot response. The aft location results in inherently stable 
operation for the bending root; vice versa for the forward location. Additional phase 
lag due to the engine position servo, sensor dynamics, and autopilot filtering (all 
desirable since phase lag implies attenuation of the high-frequency bending modes) will 
cause the direction of departure to rotate clockwise. This benefits the situation shown 


*The frequency of the TWD zero can be determined approximately by using the 
expression 
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Figure 5. Change in Root Locus Departure Angles with Sensor Location, 

First Bending Mode 

in the second diagram and deteriorates the first. Further, the modal frequency in- 
creases with flight time, meaning increasing amounts of phase lag at the bending fre- 
quency. A choice of configurations results. One can choose an aft location with mini- 
mal autopilot lag, in which case the critical flight time is toward the end of a particular 
stage of flight (the locus tends to depart upward, a marginal situation if there is some 
uncertainty in modal frequency). The other choice is a forward location with sufficient 
lag for stabilization at the beginning of a particular stage of flight (the locus departs 
more or less downward, increasing the modal stability as the flight progresses). The 
aft location implies minimal system lag and attenuation, good for the low (rigid body 
and propellant slosh) frequencies, bad for the high (higher bending) frequencies. The 
forward choice of sensor location implies the reverse. The higher-frequency modes 
may enjoy considerable stability margin and closed-loop modal damping while the 
rigid body and propellant slosh stability suffers. One may be forced into the complica- 
tion of a simultaneous switch from one sensor location and filter configuration to 
another part way through a particular stage of flight. If the frequency of this mode is 
above the TWD frequency, then the two diagrams in Figure 5 must be reversed and 
the corresponding changes made to the discussion in this paragraph. Also, it is clear 
that the 180-degree phase shift upon passing through the antinode makes a sensor loca- 
tion too close to this antinode impractical. 
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A second and more general situation is to have the rate and attitude sensors at 

different locations on the vehicle. The location of the attitude sensor is usually dic- 
tated by physical considerations, particularly if it is the inertial platform of the launch 
vehicle attitude control system. However, the rate sensor is not subject to these same 
considerations, and its compactness permits moving it to numerous other locations on 
the vehicle. 

Assuming essentially sinusoidal inputs, the rate sensor signal can be thought of 
as be mg [Kit WB ( a Rs/ CT A)J^ mes as lar 8' e as the attitude sensor signal (o^ g and o A 
are modal slopes at the rate and attitude sensor locations and tog is the frequency of 
a particular elastic mode). Thus, the influence of the attitude sensor location on the 
system stability is generally, but not always, relatively small. Under these conditions 
the remarks in this section concerning sensor locations and locus departure angles 
apply primarily to the rate sensor. Therefore, the different sensor locations could be- 
come necessary if the modal slope at the attitude sensor location was not compatible 
with elastic mode stability considerations. 

Another possibility is that of using two rate sensors, the so-called "gyro mixing" 
technique for stabilization of bending modes. This offers an additional degree of free- 
dom to the control engineer as the two rate gains can be scheduled to hold a fixed 
equivalent gain, phase, or some combination of both at bending frequencies. This 
solution might also be feasible in those situations where it is physically inconvenient to 
locate a single gyro , for example in a location that suffers from high local levels of 
vibration. And finally, it may be possible to phase stabilize two body bending modes 
at once (e. g. , 1st and 2nd body modes), which is generally beyond the capabilities of 
a single rate sensor at all flight times. Although this is not generally considered to be 
a conventional autopilot configuration, it is briefly considered in this study. 

The choice of autopilot filter is made in concert with the choice of sensor location. 
The objective is to provide proper phase shift at the frequencies of the bending modes 
being phase stabilized while maximizing attenuation at high frequencies and minimizing 
the phase lag at low (below first mode) frequencies. What constitutes proper phase 
shift at bending frequencies depends on the phase margins desired, which in turn 
ultimately depends on modal accuracy and hardware tolerance assessments. Require- 
ments for 30 degrees of phase margin and 6 db of gain margin have been used success- 
fully on Atlas -type launch vehicles. 

Generally, the stability of the propellant sloshing parasitic modes is not a major 
concern during the liftoff flight condition. Primarily, this is due to the small sloshing 
masses relative to the total vehicle mass (less than 10 percent) and the low value of 
vehicle total acceleration at this time of flight. The result is a very low level vehicle 
response to the forces and moments that are generated by the sloshing propellants. 
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Similar to the elastic modes, the propellant sloshing modes of oscillation manifest 
themselves as a series of lightly damped poles and zeros in the vehicle dynamic trans- 
fer function. Experience has shown that the coupling between these modes and the 
elastic modes is negligible if the frequency separation is a factor of 2 or greater. 
However, there is usually considerable coupling between some of the sloshing modes 
because of the standard practice of grouping, in pairs, circular cylindrical tanks of 
the same diameter. In this case the slosh pendulums will couple most strongly if they 
are located on opposite sides of the vehicle center of mass. This is the situation for 
the shuttle vehicles considered in this study. 

Another condition which is applicable to these study vehicles is the disparity be- 
tween the magnitude of the two sloshing masses. This results from the fact that the 
density of liquid hydrogen is approximately l/l6 that of liquid oxygen. Thus, the 
effect of liquid hydrogen sloshing on control system stability will be considerably less 
than the effect of liquid oxygen sloshing. 

The pole-zero sequence along the imaginary axis determines the difficulty to be 
encountered in stabilizing the slosh modes. If the sequence with increasing frequency 
is zero, pole, zero, pole, the system is relatively easy to stabilize. In the absence 
of any compensation, closing a feedback path around the transfer function results in 
loci departing directly up the imaginary axis from the poles. With the net lead from 
the combination of autopilot filtering, rate feedback, and actuator lag, the loci depart 
upward into the left-half plane. This situation is illustrated in diagrams a. andb. of 
Figure 6. 

If one of the dipoles in the pole-zero sequence is reversed the locus emanating 
from that pole will depart downward and into the right-half plane, creating a closed- 
loop sloshing instability. This situation is illustrated in diagram c. of Figure 6 . 

Such a dipole has what is termed an unstable configuration. Since such a dipole requires 
lag for stabilization (to rotate the locus into the left-half plane), while the remaining 
dipole is destabilized with increasing lag, the situation usually calls for antislosh 
baffles . 

With the pole-zero sequence of dipoles in a stable configuration, autopilot lag 
tends to destabilize the system. Since the slosh frequencies generally increase with 
increasing flight time, due to increasing axial acceleration, the root locus departure 
angles from the poles representing the slosh modes tend to rotate clockwise into the 
right -half plane; flight control system lag increases at the higher frequencies. Aggra- 
vating this situation is the decreasing vehicle mass which increases vehicle sensitivity 
to sloshing. The result is a tendency for the sloshing stability problem to be at its 
worst toward the end of a particular stage of flight. 
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Up to this point the slosh damping has been ignored. In practice it is usually very 
small, particularly if the tank has a large free surface area and is "clean" internally, 
having no internal structural members. The approximation made is quite good unless 
baffles or other internal structures are added. In this case the poles associated with 
vehicle slosh modes having considerable motion in the baffled tank tend to move to the 
left. The damping is nonlinear, increasing with increasing fluid motion. Thus, some 
slosh activity is necessary to provide the damping. If one relies on baffling to stabil- 
ize an otherwise marginally stable sloshing mode, a slosh limit cycle is the result. 

Low fluid motion (low damping) leads to an unstable root, which causes a divergent re- 
sponse at the frequency of the unstable root. As the amplitude increases , so does the 
damping and the open-loop pole and closed-loop root move toward the left. The stable 
operating point is at some limit cycle amplitude where the damping is sufficient to 
place the closed-loop root on the imaginary axis. The reverse argument shows that 
an initially damped closed-loop root will move toward the 300 axis as the modal response 
amplitude dies down, reducing the damping. Clearly the amount of baffling determines 
the limit cycle amplitude. For this study slosh damping was assumed to be zero to ob- 
tain the maximum degree of sensitivity of the control system stability to these sloshing 
modes . 

Slosh limit cycles also arise as a result of the nonlinear nature of electrohydraulic 
engine positioning servos. These arise out of the servo lag's dependence on the 
amplitude of engine motion and will result in stable or unstable limit cycle operating 
points . 

This discussion serves to indicate the consequences of using one of the following 
means for stabilizing propellant sloshing: 

a. Passive stabilization — Requires modification of the internal structure of the 
tank to increase the damping (baffles). 

b. Active stabilization — Requires modification of the flight control system via 
addition of auxiliary loops (direct slosh pressure feedback or, more commonly, 
accelerometer feedback), or network shaping. 

Either of these approaches carries its own advantages and disadvantages. Active 
compensation implies additional flight control system elements that must work to en- 
sure mission success. Further, it is not always possible to synthesize such a system. 
On the other hand, there is little if any weight penalty. Passive means are more fool- 
proof, but they carry a weight penalty. In either case, additional slosh modal damping 
mitigates inflight slosh loads that will be experienced, important since the frequencies 
are within the range of atmospheric disturbances. 
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4 DUAL FLYBACK VEHICLE ANALYSIS 


This section contains the analysis results for the dual flyback vehicle. The first 
part surveys the rigid body stability and control characteristics to establish a set of 
baseline autopilot gains for the remainder of the discussion in this section. Part two 
looks at propellant sloshing mode stability at liftoff and near liftoff flight times. The 
next two parts analyze the elastic mode stability problem for symmetric and antisym- 
metric mode boundary conditions. 

Some insight into the stability problem can be gained by examining the modal 
frequencies for the dual flyback vehicle and comparing them with frequencies of current 
operational vehicles. This frequency comparison is presented in Figure 7. 

The obviously greater modal density of the dual flyback vehicle is primarily due 
to its unsymmetrical piggy-back arrangement. The two large vehicles are elastically 
coupled at two locations along their length. With this coupling arrangement, the com- 
bined vehicle contains modes which are dominated by one vehicle or the other, as well 
as modes which reflect the interaction of the two vehicles and the elastic linkage sys- 
tem. This is in contrast to the operational vehicles, whose stages are joined in tandem. 
With these vehicles the elastic properties of the individual vehicles and the interstage 
adapters tend to weld together to form a single composite vehicle. 

A comparison between the dual flyback vehicle and the Atlas/Centaur vehicle 
illustrates this difference in the modal properties of the two types of vehicles. The 
first and second composite modes for the Atlas/Centaur vehicle exhibit the classic 
first and second body bending mode shapes. For the dual flyback vehicle the composite 
symmetric mode shapes (reference Appendix C) exhibit a booster first body bending 
mode on modes 2, 3 and 5; and modes 6, 7 and 9 contain a booster second body bending 
mode. Of course, the indicated booster mode is not dominant in all of these composite 
mode shapes, but it will be detected by the booster autopilot sensors and must be con- 
sidered accordingly in the stability analysis. 

The dual flyback vehicle mode shapes are also influenced by wing modes and 
vertical fin modes, which are not present on the other two configurations. 

There is a frequency separation of approximately a factor of two between the 
propellant sloshing modes and the first elastic modes. This fact indicates that there 
is only a small degree of coupling between these highly tuned (lightly damped) modes of 
oscillation. 
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Figure 7. Modal Frequency Comparison 

This analysis considers only autopilot sensors that are located on the booster 
vehicle. Since these sensors are required by the booster flyback capability they can 
also readily be used during the ascent flight phase. 


Rigid Body Mode 


The root locus diagrams for the rigid body modes in each of the three control 
planes are shown in Figure 8 . For these diagrams the three control planes are con- 
sidered uncoupled from each other (yaw/roll coupling effects will be discussed later). 
The damped frequency and damping ratio for autopilot gains of = l. o and K = 0. 7 
are summarized below. R 


Control Plane 


Response Parameter 

Pitch 

Yaw 

Roll 

Damped Frequency (rad/sec) 

1. 05 

1. 3 

1.6 

Damping Ratio 

0.45 

0. 45 

0. 5 
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Figure 8. Rigid Body Modes 


The damping ratios are somewhat larger than required for this time of flight, 
which indicates that the rate gain (K^) could be reduced. This would result in a de- 
crease in system phase lead (increase in lag) rotating the locus toward the right-half 
plane and decreasing the damping ratio. Lag compensation that may be required for 
elastic mode stability will produce the same type of effect. It is also apparent that the 
value of K can be less than 1. 0 in the yaw and roll control planes. 

J\ 

Propellant Sloshing Mode Stability 


Propellant sloshing stability is generally not a problem at liftoff. The sloshing 
masses are only a small percentage of the vehicle mass, and the vehicle acceleration 
is not much greater than 9. 8 m/sec 2 (1 g). For the dual flyback vehicle, the oxidizer 
(LOg) and fuel (LH^) masses are 7% and 0. 6% of the reduced vehicle mass (total mass 
less masses of slosh pendulums), and the vehicle acceleration is 13.2 m/sec 2 (1.35 g). 
These conditions are present at an approximate flight time of liftoff plus 10 sec; this time 
having been chosen to permit some buildup in sloshing masses and vehicle acceleration. 
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Figure 9 . Key Sloshing Parameters 


Figure 9 illustrates the chan- 
ges which the key sloshing param- 
eters undergo as a function of flight 
time. As long as one or both of the 

fluid levels is in the forward dome 
section of the tank (reference Ap- 
pendix C) the propellant sloshing 
natural frequencies are different. 

This frequency separation, although 
small, tends to reduce the degree of 
slosh coupling during liftoff. The 
force effectiveness and moment 
effectiveness curves indicate that 
oxidizer sloshing will have a greater 
influence on system stability than 
fuel sloshing. The different signs 
of the two moment effectiveness 
curves means that the two sloshing 
pendulum hinge points are on opposite 
sides of the vehicle center of mass. 
This condition is favorable to slosh 
coupling when the two natural fre- 
quencies are the same. 


Propellant sloshing stability at the study condition of liftoff + 10 sec is shown in 
Figure 10. These root locus diagrams are plotted in terms of damping ratio and fre- 
quency to illustrate the magnitude of tank damping that would be required to stabilize 
an unstable condition. The equations that were used to generate these loci assume 
zero natural tank damping. 

The single-order actuator model used here is adequate for the range of frequen- 
cies being considered in this analysis. The two values of K c reflect two different 
levels of effective system lag filtering that can be produced by the actuator, with 
K c = 5 depicting the higher lag condition. This varying lag is primarily the result of 
actuator nonlinearities. 


Both tanks exhibit pole above zero low residue loci. Although the hydrogen tank 
is unstable at a gain of = 1 it is a very slowly divergent condition, which requires 
less than 0. 1% damping to stabilize. This amount of damping is on the order of what is 
usually assumed for wall wiping. 

The loci in Figure 11 portray a later flight condition, liftoff + 40 sec. Now, 
both fluid levels are in the cylindrical section of each tank and the sloshing natural 
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frequencies are equal. Although the stability characteristics at a gain of = 1 are 
not appreciably different, the loci have taken a different form from that shown in 
Figure 10. The oxidizer pole has moved slightly below the zero, which indicates that 
an unstable configuration is developing. Also, this figure shows the start of extensive 
coupling between the two slosh modes. The indications are that propellant sloshing 
stability will become more critical for later flight times. 

It can be concluded that propellant sloshing is not a stability problem at liftoff. 
Therefore, for the remainder of this study, only the oxidizer slosh mode is included 
in the system equations. The reason for including this mode is to determine if it 
couples in any significant manner with the elastic modes. 


Symmetric Elastic Mode Stability 


One of the most important considerations relative to elastic mode stability is the 
location of the autopilot sensors; in particular the rate sensor. For this vehicle the 
values of the modal slopes are such that the location of the attitude sensor has only a 
secondary effect on elastic stability. Therefore, it is placed forward in the nose of the 
booster, which is the most favorable location from the standpoint of operating environ- 
ment. This analysis then examines the effects of rate sensor locations forward and aft 
of the principal first body-bending mode antinode. These sensor locations are illus- 
trated in Figure 12 in relation to the elastic vehicle model. 

For purposes of discussion in this report a distinction is made between the phrases 
"elastic mode" and "body mode". "Elastic mode" (or just "mode") refers in general to 
each of the composite vehicle modes, whereas the phrase "body mode" is used to identify 
a particular mode shape that is associated with a particular part of the total vehicle. 

For example, referring to the elastic mode shapes appearing in Appendix C (Figure 
C3) it is clear that elastic modes 2 and 3 both contain the first booster body mode shape, 
and elastic mode 9 contains the second booster body mode. Although this distinction 
is admittedly artificial, the reference to body modes more clearly describes the mode 
shape from the standpoint of the modal relationship between the autopilot sensor locations 
and the engine gimbal point. 

In the table below the symmetric vehicle elastic modes are classified according 
to the body mode shape that is exhibited on the booster vehicle. Although elastic modes 
1 and 4 exhibit elastic deflections, the modal slope has the same sign over the entire 
length of the booster. 
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BOOSTER 

STATION 



Elastic 

Mode 

Booster Body 
Mode 

Comments 

1 

(Orbiter Mode) 

Very small 
deflections 

2 

First 

Principal first 
mode 

3 

First 


4 


Large aft-end 
deflections 

5 

First 

Very small 
deflections 

6 

Second 

Predominantly a 
longitudinal mode 

7 

Second 


8 

First 

Very small 
deflections 

9 

Second 

Principal second 
mode 


For this vehicle, phase stabilizing the 
first body mode could require phase stabili- 
zation for the first five elastic modes. It is 
also possible that some of these first body 
modes can be gain stabilized, because they 
have a very small effective modal gain. 


The elastic stability problem is initially analyzed with an uncompensated autopilot 
for two rate sensor locations: one forward of the first body mode antinode at model node 
204,, and one aft of this antinode at model node 207. These locations are considered to 
be sufficiently removed from the location of the antinode to preclude a first body mode 
slope sign reversal at either location during the ascent flight phase. The root locus 
diagram for the forward rate location is shown in Figure 13, and for the aft location 
in Figure 14. 


In both of these figures the over-all stability picture is illustrated by the diagram 
on the left. The other diagrams are enlargements of specific areas of interest; namely, 
those areas which contain the loci of the individual elastic modes whose stability pro- 
perties cannot be readily deduced from the overall diagram. 
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Figure 13 shows that modes 2 and 3 must be phase stabilized since the modal gains 
are large enough to preclude gain stabilization. The other lower frequency first body 
mode (mode 5) locus also departs toward the right-half plane, but it is stable for all 
reasonable values of autopilot gain (K^). Any lag compensation which is introduced to 
stabilize modes 2 and 3 will certainly enhance the stability of mode 5. The mode 1 locus 
has a stable orientation because the modal slope at the forward rate sensor location has 
the same sign as that at the engine gimbal point. 

Mode 4 is unstable for all values of autopilot gain. It also has a low modal gain 
because of the very small modal slopes at the autopilot sensor locations. The stability 
of this particular mode is discussed in greater detail later. 

Modes 6, 7, 8 and 9 are gain stable with at least 6 dB gain margins (from = 1) 
for all phase orientations. The gain attenuation introduced by autopilot lag compensa- 
tion will increase the stability margins of these modes. 

For the aft rate sensor location (Figure 14) all mode loci exhibit stable configura- 
tions. Mode 4 is stable, despite the fact that its open-loop pole is in the right -half plane, 
because of the increased modal slope at the rate sensor location. 

It was pointed out earlier that a forward rate sensor location is most desirable 
from the standpoint of elastic mode stability toward the end of the ascent flight phase. 
With this in mind the forward rate sensor configuration is discussed in more detail. 

Figure 15 illustrates the stability conditions for a rate sensor located at model 
node 201, the location of the attitude sensor. In comparing the loci in this figure with 
those in Figure 13 the principal differences are increased modal gains and some phase 
shifts. These differences are due to the larger modal slopes at the more forward rate 
sensor location, and they are more pronounced for the high frequency modes than for 
the lower frequency modes. The changes in the second body modes are particularly 
noticeable because the rate sensor location at model node 204 is very close to a second 
body mode antinode. However, the autopilot lag compensation required to phase stabilize 
modes 2 and 3 will result in a gain attenuation factor of 3 to 4 at these second body mode 
frequencies. This comparison also illustrates the general rule that increasing gains 
have a destabilizing influence on the elastic mode loci, in that a greater reliance must 
be placed on the more uncertain technique of phase stabilization, particularly for the 
higher frequency modes. 

Returning to a rate sensor location at model node 204, Figure 16 describes the 
conditions obtained if the load torque feedback around the control actuator is removed 
(these terms are eliminated from the actuator equation). Comparing this figure with 
Figure 13 reveals that the load torque feedback effect is most pronounced at the higher 
frequencies, near the engine actuator frequency. Its dominant influence on the root 
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Figure 15 , Forward Rate Sensor (Node 201), Uncompensated Autopilot 



0.0 (X, OPEN-LOOP POLE), 
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Forward Rate Sensor (Node 204), Uncompensated 
Autopilot, Actuator Load Torque Feedback Neglected 



locus plot is to shift the location of the mode and actuator open-loop poles, in some in- 
stances toward increased stability and in others toward decreased stability. In particu- 
lar, note the destabilizing effect that the load torque feedback has on mode 4. 

This figure illustrates that the stability of mode 4 will be very much influenced 

by the effective spring rates for the actuator backup structure for each of the twelve 
engines and the friction terms and nonlinearities inherent in the engine actuator sys- 
tems. Since none of these effects was included in this simplified linear analysis, this 
discussion only indicates a potential problem area that must be carefully considered in 
the design phase. Generally, the stability of this type of mode can be modified to some 
degree by the control the designer has over actuator system requirements and actuator 
backup structure spring rates. 

If stability compensation is added to the autopilot for phase stabilizing modes 2 
and 3 (model node 204 rate sensor location), the results are shown in Figure 17. The 
compensation consists of a forward loop quadratic lag filter with a natural frequency of 
15 rad/sec and critical damping (£ = 0.707). Now all of the elastic mode loci (again, 
with the exception of mode 4) have stable closed-loop roots for the range of autopilot 
gains of interest. 

Although mode 1 has adequate stability margins, the introduction of the filter has 
had a destabilizing influence on it. This mode is predominantly an orb iter mode, whose 
frequency changes very little with flight time. Therefore, as the frequencies of the 
other modes increase with flight time, the amount of autopilot lag filtering can be re- 
duced with a resulting increase in the stability of this mode. 

For this study it was practical to consider only briefly modal dispersions on 
frequency. The root locus plots for plus and minus ten per cent variations in elastic 
mode frequency are shown in Figures 18 and 19 using the compensated autopilot con- 
figuration. With some small adjustments to the filter natural frequency, adequate 
stability margins can be attained for all elastic modes, excluding mode 4. Certainly 
the conventional autopilot techniques are adequate to cope with these frequency 
dispersions. 

If the stability of a mode with characteristics similar to those of mode 4 cannot 
be attained by careful actuator system design, a multiple rate sensor autopilot con- 
figuration may be required. This concept was very briefly examined by considering 
an autopilot with two rate sensors, one located forward of the first body mode antinode 
and one mounted aft of this antinode (model nodes 204 and 207). The influence of each 
sensor can be controlled by adjusting its rate gain. For this example the forward sen- 
sor has a rate gain (K^) of 0. 2 and the aft sensor a gain of 0. 5. Thus, the effective 
rigid body gain is still 0. 7. The root locus plot for this configuration is shown in Figure 
20- Now, all elastic modes exhibit stable characteristics, including mode 4. 
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Figure 17. Forward Bate Sensor (Node 204), Compensated Autopilot 
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Figure 18. Forward Rate Sensor (Node 204), Compensated Autopilot, Elastic 
Mode Natural Frequencies Increased 10 Per Cent 
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Figure 19. Forward Rate Sensor (Node 204), Compensated Autopilot, Elastic 
Mode Natural Frequencies Decreased 10 Per Cent 




In conclusions this analysis has demonstrated the feasibility of using conventional 
autopilot techniques to provide the pitch control plane of the dual flyback vehicle with 
adequate stability and performance characteristics at the liftoff flight condition. 


Antisymmetric Elastic Mode Stability 


Antisymmetric elastic mode stability is analyzed in two phases. In the first phase 
a planar analysis is conducted for the yaw and the roll control planes separately. This 
initial planar approach permits greater visibility of the principal yaw and roll body 
mode stability properties by eliminating the effects of the vehicle dynamic model coupl- 
ing terms. Then, in the second phase the coupled system is examined in terms of the 
results of these planar analyses. 

In the table below the antisymmetric vehicle elastic modes are classified accord- 
ing to the body mode shape exhibited on the booster vehicle. The yaw and roll com- 
ponents are listed separately. 


Elastic 

Booster Body Mode 


Mode 

Yaw 

Boll 

Comments 

1 

First 


Predominantly a wing mode 

2 

First 


"Scissor mode" between booster and orbiter 

3 

First 


Principal first yaw mode 

4 

First 



5 

First 



6 

Second 

First 

Principal second yaw mode 

7 

Second 



8 

Second 



9 


First 

Principal first roll mode 


Unlike the symmetric modes, this set of modes does not contain a mode that is 
almost entirely composed of booster aft end deflection. This could be due to the differ- 
ent total vehicle elastic properties between the longitudinal and lateral planes. Another 
factor to consider is the difference in elastic models. The antisymmetric modes are 
obtained from a three-dimensional representation of the booster vehicle, whereas the 
symmetric modes use only a flexible beam model. This more detailed antisymmetric 
model may have a distributive effect on a "localized" mode shape. 

Unless noted, the mode shape data (deflections and slopes) are expressed in terms 
of equivalent booster centerline parameters. 
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The rate and attitude sensor locations 
examined here are illustrated in Figure 21 
in relation to the elastic vehicle model. 

Yaw Planar Analysis . - The root locus 
diagrams for rate sensor locations forward 
and aft of the first body mode antinode are 
shown in Figures 22 and 23. The forward 
location is at model node 740 and the aft 
location is at model node 1040 (the gimbal 
point). The only attitude sensor location 
considered is forward in the booster nose 
(model node 600). 

Figure 22 shows that modes 2 and 3 
are unstable (closed-loop roots in the 
right-half plane). Also, the modal gain 
is large enough that these two modes 
cannot be gain stabilized. Mode 1 is 
represented by a stable dipole configura- 
tion with a very small residue. In this 
configuration Mode 1 has a negligible 
influence on the stability properties of 
the other system modes. 


1000 

1020 

1040 


Although the loci for modes 4, 5 and 6 
depart toward the right-half plane, the 
closed-loop roots for = 1 are in the left- 
half plane indicating stability. Also, modes 
7, 8 and 9 are completely gain stable (i. e. , 
these modes have at least a 6 dB gain 
margin for all phase orientations). Any 
lag compensation introduced to stabilize 
modes 2 and 3 will certainly increase the 
stability margins of modes 4 through 9. In this figure the second body mode loci 
depart toward the right-half plane. This condition is primarily due to the location 
of the rate sensor just aft of the second mode forward antinode. A more forward rate 
sensor location will result in these loci departing toward the left. 


Figure 21. 


Autopilot Rate and Attitude 
Sensor Locations 


For the aft rate sensor location (Figure 23), all elastic mode loci exhibit stable 
configurations. 
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Figure 22. Planar Yaw Control, Forward Rate Sensor (Node 740) 
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Figure 23. Planar Yaw Control, Aft Rate Sensor (Node 1040) 



Returning to the forward rate sensor location (model node 740) Figure 24 illus- 
trates the effect of adding autopilot compensation to provide phase stabilization for 
modes 2 and 3„ The compensation consists of a forward loop quadratic lag filter with 
a natural frequency of 10 rad/sec and critical damping. Now, all of the elastic mode 
loci have stable closed-loop roots for a gain value of = 1. Also, modes 5 through 9 
are completely gain stable. 

Although the introduction of this autopilot filter has a destabilizing influence on 
the propellant sloshing and rigid body modes, it does not produce unacceptable results. 
The propellant sloshing locus still retains its stable configuration (see Figure 22). At 
K^ = l. 0 the rigid body damping ratio has been reduced to 0. 35 from 0. 45, and the 
damped frequency has increased to 1. 5 rad/sec from 1. 3 rad/sec. 

It is shown in the coupled system analysis that the coupling effects do not appre- 
ciably change the results obtained from this planar analysis. The differences that are 
noted in the coupled root loci are small from the standpoint of overall vehicle stability. 
Therefore, the effect of elastic mode frequency dispersions can be analyzed in this 
planar condition and the results, if significant, can be interpreted in terms of the 
coupled system. 

Using the compensated forward rate sensor configuration, the root locus diagrams 
for a plus ten percent and a minus ten percent change in elastic mode natural frequencies 
are given in Figures 25 and 26 . These figures illustrate that all loci still have stable 
closed -loop roots at a gain of K = 1. 0. The somewhat marginal stability qualities 

jlx 

exhibited by mode 2 for the decreased frequency conditions are easily improved by 
small adjustments to the autopilot compensation. 

Roll Planar Analysis . - The root locus diagrams for forward and aft roll rate 
sensor locations are presented in Figures 27 and 28. The forward rate sensor location 
is in the booster nose at model node 600, and the aft location is at model node 1040 
(the gimbal point). For both rate sensor locations, the roll attitude sensor is positioned 
in the booster nose. 

In comparing these two figures, only modes 6 and 9 exhibit the approximately 
180-degree phase shift associated with first body modes when the rate sensors are on 
opposite sides of the antinode. Although the other modes exhibit roll deflections, 
there is no sign change between the aft and forward sections of the vehicle. Hence, the 
loci representing these modes assume the inherently stable configuration characteristic 
of a first mode with an aft-mounted rate sensor. 

For both rate sensor locations all elastic mode loci have stable closed-loop roots 
at Ka = 1. 0. Also, in both figures the mode 9 locus exhibits the smallest stability 
margins. Although this mode could be classified as a higher frequency mode, it 
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Figure 24. Planar Yaw Control, Forward Rate Sensor (Node 740), 
Compensated Autopilot 
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Figure 25. Planar Yaw Control, Forward Rate Sensor (Node 740), 
Compensated Autopilot, Wbejjt) + 10% 
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Figure 26. Planar Yaw Control, Forward Rate Sensor (Node 740), 
Compensated Autopilot, - 10% 
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Figure 27. Planar Roll Control, Forward Rate Sensor (Node 600) 
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Figure 28. Planar Roll Control, Aft Rate Sensor (Node 1040) 



represents the principal first body mode and has sufficiently high modal gains to re- 
quire phase stabilization at this early flight time. The preferred orientation of this 
locus is attained with the forward rate sensor location because of the tendency of the 
locus to rotate clockwise with increasing modal frequency. 

With this large complex and unsymmetrieal structure, it might be anticipated 
that greater consideration should be given to local structural deformations in deter- 
mining satisfactory autopilot sensor locations for this vehicle than has been the 
practice on other vehicles. It then follows that sufficiently detailed elastic models 
must be required to evaluate these effects. The use in this part of the study of a 
partial three-dimensional elastic model permits a brief evaluation of this situation, 
particularly as it pertains to the model. 

For the mode shapes that were derived from this three-dimensional model, the 
roll characteristics have the largest differences between the equivalent center- 
line representation and the representation of each node around the periphery of the 
vehicle. These differences encompass both the magnitudes and the signs of the roll 
deflections . 

The root locus diagram for the condition where the roll rate sensor is located 
at one of these peripheral nodes is shown in Figure 29. The particular node repre- 
sented in this figure is on the underside of the vehicle, just aft of the canard at S. N. 
2170 (model node 780-7). The roll modal data reveals that Figure 28 can be used 
for a comparison between this peripheral rate sensor location and a location that 
reflects an equivalent centerline position, particularly relative to loci departure angles. 

The loci departure angles for several elastic modes have undergone an approxi- 
mately 180-degree phase shift at the underside rate sensor location, and modes 2, 4, 

6, 7 and 9 require some degree of phase stabilization. In particular, mode 4 at 11. 7 
rad/sec requires greater than 90 degrees of phase lead or phase lag. From the stand- 
point of conventional autopilot techniques, this is not a practical location for the roll 
rate sensor. 

A review of the roll modal data at the model node located on the topside of the 
vehicle (at model node 780) reveals that these data are very similar to the equivalent 
centerline data. Thus, a rate sensor located at this topside node will have similar 
stability characteristics to those shown in Figure 28. 

The effects of using the aft rate sensor configuration with dispersions in the 
elastic mode natural frequencies of plus and minus ten percent are shown in Figures 
30 and 31. These figures illustrate that all loci still have stable closed-loop roots at 
a gain of -1.0. 


50 


LEGEND 


K - 0.0 (x, OPEN- LOOP POLE), 

0.2, 0.6, 1.0 (*), 1.4, 1.8, 
2.0, 3.0 AND «> (0, ZERO) 

ACTUATOR POLE 




Figure 29. Planar Boll Control, Underside Rate Sensor (Node 780-7) 
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Figure 33. Planar Roll Control, Aft Rate Sensor (Node 1040), to 


Yaw/Boll Coupled Analysis . — With a multiloop multivariable control system, 
each response can be related to each input by a transfer function. The common de- 
nominator for all these transfer functions is the system characteristic equation. As 
with single variable systems, stability information can be obtained from the root loca- 
tions of this equation. 

A limitation of the root locus method is that only one varying parameter can be 
considered per locus. Therefore, to analyze the stability of this coupled yaw/roll 
system, the forward loop gain, K^, in one of the control planes is fixed, and the locus 
is generated by varying the corresponding gain in the other control plane. For example, 
if Ka in the roll autopilot is fixed and in the yaw autopilot is varied, then the re- 
sulting loci will describe the stability characteristics of the yaw control system in 
the presence of the roll-plane coupling effects. The two fixed values of gain used here 
are 0. 0 (open loop) and 1. 0 (assumed nominal operating point). The fixed value of 0. 0 
permits an evaluation of the coupling effects that are not due to the other control plane 
autopilot (i. e. , intermodal coupling and coupling through the engine gimbal actuators). 

For this analysis, the forward rate sensor configuration with autopilot stability 
compensation (lag filter with natural frequency of 10 rad/sec and critical damping) is 
used for yaw control, and the aft rate sensor configuration is used for roll control. 

Since the forward roll rate sensor configuration demonstrated satisfactory stability 
characteristics, it could be used in place of the aft rate sensor. 

Figure 32 shows the stability characteristics for the yaw control system with 
the roll autopilot gain set equal to zero. The vehicle dynamic and load torque feed- 
back coupling effects can be evaluated by comparing this figure with Figure 24. The 
rigid body locus is unchanged, and the propellant sloshing locus has retained its stable 
orientation. The mode 1 locus shows a large effective gain increase (compare with 
Figure 22), but the K^ = l point still has a greater than 6 dB gain margin. This wing 
mode appears to be interacting strongly in the region of the engine thrust structure. 

The other lower frequency mode loci exhibit only very small differences in the range 
of gain values that are of interest here. 

The higher frequency modes in the vicinity of the tail-wags-dog zero exhibit 
significant differences, as indicated by the changes in the zero locations. However, 
these modes are still completely gain stable. Since the location of the TWD zero is 
now being influenced by both the yaw and roll actuator systems, it also has shifted 
relative to its planar location. 

The stability characteristics for the roll control system with the yaw autopilot 
gain set to zero are illustrated in Figure 33. In a comparison with Figure 28 only very 
small differences are noted. 
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Figure 33. Coupled Yaw/Roll, Vary Roll Autopilot Gain, Yaw Gain = 0. 0 


The completely coupled system is illustrated in Figure 34 with a variable yaw 

autopilot forward loop gain and a fixed value of 1. 0 (nominal operating point) for the 
roll autopilot gain. The loci in this figure have the general shape shown in Figure 32, 
and the open loop poles (the Xs) for the bending mode loci are located at the position 
of the closed-loop roots (K A = 1. 0 points) in Figure 33. Figure 34 shows that the loci 
for Modes 1, 4, 5, 6, 7, 8 and 9 are gain stable with at least 6 dB gain margin for all 
values of phase, and that Modes 2 and 3 are phase stable with phase margins greater 
than or equal to thirty degrees. 

Figure 35 is the corresponding plot for a variable roll autopilot forward loop gain 
and a fixed value of 1. 0 (nominal operating point) for the yaw autopilot gain. Here, the 
loci have the general shape shown in Figure 33, and the open-loop poles for the bending 
mode loci are located at the position of the closed-loop roots in Figure 32. All of the 
bending mode loci in Figure 35 are adequately stabilized with Modes 3, 6, and 8 gain 
stable for all phase orientations. 

In conclusion, this analysis has demonstrated the feasibility of using conventional 
autopilot techniques to provide lateral (yaw/ roll) control for the dual flyback vehicle 
with adequate stability and performance characteristics at the liftoff flight condition. 
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Figure 34. Coupled Yaw/Roll, Vary Yaw Autopilot Gain, Roll Gain - 1. 0 





5 SOLID MOTOR VEHICLE ANALYSIS 


This section contains the analysis results for the solid motor vehicle. The first 
part surveys the rigid body stability and control characteristics to establish a set of 
baseline autopilot gains for reference in the remainder of this section. This part 
also discusses propellant sloshing mode stability. The last two parts analyze the 
elastic mode stability problem for s y mmetric and antisymmetric mode boundary 
conditions. 


Figure 36 compares the solid motor vehicle modal frequencies with those for 
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Figure 36. Modal Frequency Comparison 

the dual flyback vehicle. The solid motor vehicle has higher first elastic mode fre- 
quencies and the resulting frequency separation of greater than a factor of two be- 
tween the propellant sloshing modes and the lower elastic modes greatly minimizes 
the coupling potential between these modes. 

The decrease in the tail-wags -dog frequency indicates that this singularity will 
have a greater influence in determining the stability characteristics of the solid 
motor vehicle first body modes than it did for the dual flyback vehicle. 
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Control forces for this vehicle can be obtained from three sources: the orbiter 
liquid propellant engines, the solid rocket motors, and the orbiter aerodynamic sur- 
faces. Current shuttle concepts favor using thrust vector control on the orbiter en- 
gines with the possibility of using the orbiter aerodynamic surfaces if additional con- 
trol forces are required during the high dynamic pressure region of flight. Thrust 
vector control on the solid motors is not considered a viable alternative at this time 
because of the cost and weight penalties that are associated with its implementation. 

The baseline configuration for this study uses only thrust vector control of the 
three orbiter engines. Aerodynamic surfaces are not considered because of their 
relative ineffectiveness at the liftoff flight condition. Thrust vector control on the 
solid motors is briefly considered to demonstrate its effect on elastic mode stability. 


Rigid Body and Propellant Sloshing Modes 

There are two areas of concern relative to the magnitude of the control forces 
available on the solid motor vehicle during the atmospheric portion of flight. The 
first involves the control authority (gimbal angle capability) required to satisfactorily 
maintain the vehicle's desired heading. The second pertains directly to this study 
and involves the stability and performance characteristics of the vehicle control 
mode. 


Some insight into the second area at the liftoff flight condition can be obtained 
by examining the control moment effectiveness ratios (control thrust times control 
moment arm divided by moment of inertia) for each of the control planes; 

Solid Motor Vehicle 


Control 

Plane 

Orbiter 
TVC Only 

Orbiter TVC + 
Solid Motor TVC 

Dual Flyback 
Vehicle 

Pitch 

0.31 

1.6 

1.8 

Yaw 

0.27 

1.4 

1.8 

Roll 

0.11 

2.4 

3.0 


If the orbiter engines and the solid motors are both used for control, the solid 
motor vehicle ratios are comparable to those for the dual flyback vehicle. However, 
if only the orbiter engines are considered, the ratios are considerably less than 
those for the other two configurations. This implies that the "baseline" solid motor 


62 


vehicle (orbiter TVC only) will require higher autopilot gains to obtain satisfactory 
control mode response characteristics . 

The control mode gains can be determined from the pitch, yaw and roll control 
mode root loci shown in Figure 37. For a minimum response with a damping ratio 
not less than 0.3 and a damped natural frequency of approximately 1 rad/sec, the 
autopilot forward loop gains are: pitch = 4. 0, yaw = 4. 0, and roll = 10. 0. The auto- 
pilot rate gain (Kr) for each control channel is 0. 7. 

Propellant sloshing will have a negligible effect on the overall vehicle stability 
characteristics at the liftoff flight condition. The oxidizer and fuel sloshing masses 
are only 0. 9% and 0. 3% of the vehicle mass. The tank levels that correspond to this 
time of flight are shown in Figure D1 in Appendix D. It has already been shown in 
Figure 36 that the lowest elastic mode frequencies are approximately a factor of three 
removed from the sloshing natural frequencies, indicating a negligible coupling poten- 
tial between these very lightly damped modes. 

The sloshing natural frequencies are illustrated in Figure 38 as a function of 
flight time. The two sloshing modes have different natural frequencies because the 
oxidizer fluid level is in the conical section of the tank during these early flight times. 
This frequency separation and the relatively small fuel sloshing mass will tend to fur- 
ther minimize the coupling potential between these sloshing modes. 

Based upon these facts and the results of the dual flyback vehicle analyses it is 
assumed that propellant sloshing effects can be neglected at this time of flight. In 
order to evaluate this assumption and detect any potential surprises the vehicle 
model used in this section includes the oxidizer sloshing mode. 

Symmetric Elastic Mode Stability 

From the standpoint of stability and control, the solid motor space shuttle con- 
figuration represents a major design departure from present operational vehicles, 
and even from the dual flyback vehicle. Four large bodies, three of which have quite 
different structural properties are coupled together in a parallel arrangement. This 
produces a situation where there can be a large amount of interaction between the in- 
dividual elements of the vehicle through the Connecting linkage systems. 

Some indication of this interaction can be obtained by examining the elastic 
mode shapes contained in Appendix D (Figure D3), The body bending mode shapes 
are summarized in the following table. 
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Figure 37. Rigid Body Root Loci 



Elastic 


Body Mode 



Mode 

Orbiter 

Tank 

SRM 

Comments 

1 



"Bridge" 

Orbiter pitching and plunging on 
links 

2 




Orbiter pitching motion on links 

3 

First 


First 

Dominant orbiter mode 

4 

First 

First 

First 

Orbiter opposing tank and 
dominant SRM mode 

5 




Large orbiter forward end 
deflections 

6 

Second 



Dominant ox’biter forward end 

7 

Second 

First 

Second 

Dominant SRM mode 

8 

Second 

First 


Dominant tank mode 


Only those body modes that exhibit all of the classical properties of the indicated 
mode shape are listed. In the first elastic mode the SRM element mode shape rep- 
resents the classic first mode for an elastic body that is pinned at both ends. 

There is an important difference between the orbiter mode shape in elastic mode 
5 and the orbiter first body mode shapes in elastic modes 3 and 4. In mode 5 the sign 
relationship between the modal slope along the length of the element and the modal 
deflection at the engine gimbal point is the opposite of the corresponding relationship 
for a first body mode. If the autopilot sensors are mounted on the orbiter, this sign 
difference will appear in a root locus diagram as a 180 degree phase shift between the 
locus departure angles from the mode 4 and 5 open loop poles. Considering the small 
frequency separation between these two modes , this condition will complicate the 
elastic mode stability problem, particularly if one or both of these modes is dependent 
on some form of phase stabilization. 

Although the preferred location for the autopilot sensors is on the orbiter , the 
stability complications that developed indicated that it might be more advantageous 
to locate these sensors on the tank or the SRM. The sensor locations examined in 
this analysis are indicated in Figure 39. 

The root locus diagrams for orbiter rate sensor locations forward and aft of the 
first body mode antinode are shown in Figures 40 and 41. The forward location is at 
model node 4 and the aft location is at model node 7. For both cases the attitude sensor 
is located in the orbiter nose (model node 1). 
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Figure 38. Solid Motor Vehicle Propellant Sloshing Natural Frequencies 
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Figure 39. Autopilot Sensor Locations 


The stability consequences of the high autopilot forward loop gain are clearly 
evident in Figure 40. The high modal gain forces the closed-loop roots to be signifi- 
cantly displaced from the open-loop poles. This makes modal phasing particularly 
critical, even for the higher frequency modes. Theoretically, the modal gain attenu- 
ation required to stabilize this configuration (approximately -17 dB) could be obtained by 
locating the rate sensor very close to the first body mode antinode . However , the 
uncertainties inherent in the location of an antinode for all flight times and the exist- 
ence of two first body mode antinodes (one for each first mode) eliminates this solution 
from practical consideration. 

The location of the tail -wags -dog zero at approximately 19 rad/sec has placed 
this singularity below the orbiter first body mode natural frequencies , but above the 
two lowest frequency elastic modes. As a result of this location the loci emanating 
from the first body mode open loop poles are departing upward, instead of the down- 
ward orientation that is customarily associated with a forward rate sensor configuration. 

The highly unstable condition at the nominal operating gain precludes the use of 
conventional autopilot techniques to stabilize this configuration. Adequate gain 
attenuation cannot be practically obtained and phase compensation is limited by the 
opposing nature of the mode 4 and mode 5 loci and the small frequency separation 
between these modes. 

For the aft rate sensor location (Figure 41) the effect of the high autopilot for- 
ward loop gain is evidenced by the nearness of the nominal operating gain to the 
terminal zero. However, except for the locus from the mode 3 open-loop pole, all 
elastic mode loci are completely inside the left-half plane . 

Although modes 2 and 3 have similar modal parameters relative to the rate 
sensor location and the engine gimbal point, the presence of the TWD zero has re- 
sulted in an approximately 180 degree difference in their locus departure angles, 
and an unstable system. 

The small frequency separation between these two modes makes phase stabilization 
extremely difficult. An example of this is shown in Figure 42 where a single lag filter 
with a 0. 008 second time constant has been added to the autopilot forward loop. 

Although this filter contributes only ten degrees of phase lag at 24 rad/sec it is 
sufficient to cause the loci to "swap" zeros. Now, any additional phase lag will in- 
crease the system instability by producing a clockwise rotation of the mode 2 locus. 

The effects of modal frequency uncertainties in this unstable region are simulated 
by varying the natural frequency of mode 2 by minus 2 rad/sec and by plus 2. 9 rad/sec. 
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Figure 40 . Forward Rate Sensor on Orbiter (Node 4) 
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Figure 41. Aft Rate Sensor on Orbiter (Node 7) 
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Figure 42. Aft Rate Sensor on Orbiter (Node 7) , Compensated Autopilot 

The root locus diagrams for these two conditions, along with the one for the "nominal" 
value of frequency, are shown in Figure 43. 

Increasing the frequency separation between the mode 2 open-loop pole and the 
TWD zero increases the system stability. In both cases the system is now more 
amenable to conventional autopilot stability compensation, but the high gain still makes 
the attainable degree of stability marginal . 

Thus far conventional autopilot techniques have been considered for a configuration 
which has the autopilot sensors mounted on the orbiter and uses only orbiter engine 
TVC. The high modal gains and the closeness of the elastic mode frequencies limits 
the ability of these techniques to provide adequate stability compensation for tins con- 
figuration. 

Although the two rate sensor concept is not, strictly speaking, a conventional 
autopilot technique , the extra degree of freedom that it provides the controls designer 
could be very useful for this vehicle. Figure 44 shows the results for a two rate 
sensor configuration where the rate gain for the forward rate sensor at node 4 is 0.2 
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Figure 43. Aft Rate Sensor on Orbiter (Node 7), Vary Mode 2 Natural Frequency 
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Figure 44. Two Rate Sensors on Orbiter, Forward at Node 4 and Aft at Node 7 





and 0,5 for the aft rate sensor at node 7. The autopilot attitude sensor is located in 
the orbiter nose at node 1 . 

The rate gains are weighted in favor of the aft rate sensor because of the highly 
unstable condition produced by the forward rate sensor. The result is that the loci near 
the TWD zero are oriented in much the same way as they are for the aft rate sensor 
configuration (see Figure 41). However, the two rate sensor configuration produces 
a noticeable stability improvement by introducing some modal gain attentuation into 
the mode 2 and mode 3 loci. A small amount of lag filtering can now be used to shift 
the mode 3 closed-loop root into the left-half plane. 

Although system stability can be attained with this autopilot configuration , the 
gain and phase margins are not adequate, particularly for the modes near the TWD 
zero. The root loci orientations in this region are sensitive to the modal parameters 
and large stability margins are required to compensate for modal uncertainties. 

To round-out this symmetric mode analysis two other control configurations will 
be briefly analyzed. The first will consider mounting the autopilot sensors, principally 
the rate sensor, on the other vehicle elements (the external tank and the SRM). The 
second configuration will add SRM TVC to the orbiter TVC. 

Figure 45 (a) shows the root locus diagram for a rate sensor mounted on the 
tank near the orbiter forward linkage system (node 12) and an attitude sensor located 
in the orbiter nose (node 1). All modes exhibit stable characteristics with the exception 
of mode 3 (orbiter first body mode) and mode 8 (tank first body mode). 

When the two autopilot sensors are located on different vehicle components it is 
possible that the attitude signal may dominate the rate signal. It is easy to verify 
from the mode shape in Figure D3of Appendix D that this is the situation for mode 3. 
This also suggests that this mode could be stabilized by filtering only the attitude 
signal, particularly if the stability of the neighboring modes is not dependent on this 
signal. 

Figure 45 (b) illustrates the effect of placing a critically damped quadratic lag 
filter, with a natural frequency of 10 rad/sec, in the attitude signal path. Mode 3 has 
been stabilized without destabilizing any of the surrounding modes, although the 
orientation of some of them (e.g. , mode 6)has been significantly altered. This leaves 
only the first tank mode in an unstable configuration. Because of the high modal gains 
this mode must be phase stabilized by introducing approximately ninety degrees of 
phase lag into the autopilot forward loop. This large amount of phase lag will make 
the modes at 30 and 25 rad/sec marginally stable. 
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The high signal attenuation factor introduced by the quadratic lag filter in the 
attitude sensor feedback path minimizes the influence of this sensor on the stability 
of modes above 20 rad/sec. This condition produces approximately the same results 
obtained if the attitude sensor is placed at the same location as the rate sensor. 

Although the rate sensor position in Figure 45 is near the middle of the tank 
element, it is aft of the tank first body mode antinode. A sensor position near 
the orbiter aft linkage system (node 15) was also analyzed with results very similar 
to those for the mid-tank position. The only difference of any consequence to the 
system stability is the mode 8 gain increase for the aft sensor location. This modal 
gain increase increases the phase compensation required for stability, making this 
aft sensor location less desirable than the mid-tank position. 

Locating the rate and attitude sensors on the tank near the SRM forward linkage 
system (node 10) produces the results shown in Figure 46. Now the rate sensor is 
forward of the first tank mode antinode and mode 8 departs into the left -half plane. 

All loci exhibit stable configurations except mode 4. Unfortunately, modes 3 and 4 are 
in a basically unstable configuration. Introducing the lag compensation required to 
stabilize mode 4 will cause the two loci to "swap" zeros, with the result that the 
system will still be unstable. 

Locating the rate sensor on the forward end of the SRM (node 17) and the attitude 
sensor in the orbiter nose produces the root loci configuration shown in Figure 47. 
Modes 1 , 3 and 7 exhibit unstable orientations and filtering the attitude sensor signal 
does not alter this situation. All three modes must be phase stabilized, but the con- 
flicting phase stabilization requirements between adjacent modes rules out this 
approach. 

The root locus diagram for the rate sensor on the aft end of the SRM (node 22) 
and the attitude sensor in the orbiter nose is shown in Figure 48(a). Modes 2, 3, 4 and 
7 exhibit unstable orientations . If the attitude sensor is moved to the rate sensor 
location on the SRM (Figure 48(b)) only modes 2 and 7 are still unstable. However, any 
phase compensation introduced to stabilize these modes will destabilize either mode 
3 or mode 4. 

One of the major difficulties in attempting to solve the elastic mode stability 
problem for this vehicle is the destabilizing influence of high modal gains. These 
high gains, for the most part, are a result of the relatively high autopilot forward 
loop gain required to achieve satisfactory control-mode stability and response char- 
acteristics. To decrease this gain, it is necessary to increase the magnitude of the 
control forces. 
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Figure 47. Rate Sensor Forward on SRM (Node 17), Attitude Sensor in Orbiter Nose (Node 1) 
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To a first-order approximation, the effect of a control force increase on elastic 

mode stability generally can be evaluated with the nominal control thrust root locus 
diagrams. Simply shift the closed-loop root along the same elastic-mode locus to the 
new value of autopilot forward-loop gain (K^) determined from the control mode re- 
sponse requirements. This approximate technique is particularly applicable when 
already active controllers are being considered; such as an increase in the thrust 
of the orbiter engines. 

For this vehicle increased control thrust can be obtained by adding a TVC system 
to the solid rocket motors. Since this concept involves the addition of a different ac- 
tuator system and new force application point on the vehicle, the SRM TVC system is 
incorporated into the dynamic equations. A secondary injection system is modeled as 
a single-order lag with a time constant of 0. 062 sec. For the purposes of this study, 
the vehicle mass properties are not adjusted to reflect the added system, and the engine 
command signal applied to the orbiter gimbal system is also used for the SRM TVC 
system. 

Figure 49 shows the root locus diagram for autopilot sensors located on the orbiter 
(aft rate sensor at model node 7 and forward attitude sensor at model node 1). The en- 
largement of the control mode locus reveals that the minimum response characteristics 
can be obtained with an autopilot forward loop gain (K^) of 0. 6, as opposed to the gain 
of 4. 0 without SRM TVC. All of the elastic mode closed-loop roots (K^ =0.6 points) 
have assumed stable orientations with the exception of the slightly unstable mode 2. 

This instability reflects the complications that are introducted by the TWD zero from 
the orbiter engine gimbal system. These complications will be present for all practical 
values of gain. 

A comparison between Figure 49 and Figure 41 reveals that the reduced autopilot 
gain has improved the overall system stability and enhanced the potential for being able 
to provide adequate stability compensation. Of course, it has not eliminated all of the 
stability problems associated with this vehicle. 

This comparison also illustrates that significant changes occur in some locus 
orientations. Thus, the approximate method of estimating gain reduction effects can 
only be used to obtain a rough first-order approximation. Even with this limitation, 
it can be seen in Figure 40 that a reduction in gain to = 0. 6 will result in a signif- 
icant improvement in the system stability of this forward rate sensor configuration. 

This symmetric elastic mode stability analysis serves to delineate the diffi- 
culties associated with using conventional autopilot techniques to provide system sta- 
bility for this study vehicle configuration. The feasibility of using these techniques 
for this vehicle is not resolved by this analysis. The best results obtained with 
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a single rate sensor configuration were very marginal. Although the two-rate sensor 
configuration could provide positive system stability, it did not provide really satis- 
factory margins of stability. 


Antisymmetric Elastic Mode Stability 


As described in Section 3, the antisymmetric elastic mode stability is analyzed 
in two phases. The first phase consists of a planar analysis of both the roll and the 
yaw control planes separately. This planar analysis is artificial in the sense that the 
modal data consists of the yaw and roll components from the coupled antisymmetric 
modal data. Then the second phase analyzes the complete coupled system. 


In the table below the antisymmetric elastic modes are described in relation to 
the body modes that appear on each of the three different vehicle elements. The yaw 
and roll components of each elastic mode are listed separately. 


Lateral (Yaw) 
Projection Body Mode 
Orbiter Tank SRM 


Roll Projection 
Body Mode 

Orbiter Tank SRM 


First First 
First First 


Elastic Mode 


1 

2 

3 


4 


First 

5 


First 

6 


First 

7 


First 

8 


First 

9 

First 

First 

10 

First 

First 

11 

First 

First 

12 




First Second 
First 

Second Second 
Second 
Second 
Second 

Second 


Comment 


Dominant coupled tank/ 
SRM first yaw mode. 
Dominant orbiter first 
roll mode. 

Dominant tank first 
yaw mode. 

Dominant orbiter 
roll mode. 

Dominant orbiter first 
yaw mode. 


First 

First 

First 

First 

First 

First 

First 

First 


For the first seven elastic modes the orbiter lateral bending appears essentially 
as a rigid body interacting with the other vehicle elements through its interconnect- 
ing linkage system. Finally in mode 8 some significant lateral elastic motion appears 
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that is in opposition to the dominant tank first-yaw mode. The shape and construction 

of the orbiter vehicle produces a much stiffer structure in yaw than in pitch. This is 
reflected in the high natural frequency of the orbiter first yaw body mode at 7.4 Hz, as 
compared to the first pitch body mode frequency at 3. 8 Hz. 


The autopilot rate and attitude sensor locations used in this study are shown in 
Figure 50 relative to the vehicle elastic model. 


LENGTH STATION 
(M) NUMBER 


YAS - YAW ATTITUDE SENSOR 
RAS - ROLL ATTITUDE SENSOR 



GIMBAL 

POINT 


COMPONENTS 


Figure 50. Autopilot Sensor Locations 
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Yaw Planar Analysis , - The root locus diagrams for rate sensor locations for- 
ward and aft of the first orbiter body mode antinode are shown in Figures 51 and 52. 

The forward location is at model node 14 and the aft location is at model node 6. The 
only attitude sensor location considered at this time is forward in the orbiter nose 
(model node 22). Several alternate locations, both on the orbiter and on the tank, were 
examined briefly. 

Figure 51 shows that mode 10 is unstable. This is the dominant orbiter first yaw 
mode, and the modal gain is large enough to preclude gain stabilization. Several modes 
(1,2, 7, 8, 9, 11 and 12) have very small residues and are completely gain stable. Modes 
3,4,5 and 6 depart somewhat into the left-half plane and are stable for all gains. Filter- 
ing to affect mode 10 (at 45 rad/sec) could cause problems with mode 5 (at 21.2 rad/sec) 
and mode 6 (at 26 rad/sec). 

Mode 10, with its open-loop pole at 45 rad/sec poses a severe problem because 
of the proximity of the engine actuator pole at 58. 5 rad/sec. Any filtering to attempt 
to rotate the locus of mode 10 counter clockwise will result in "loci swapping"; i. e. , 
the actuator locus becoming unstable and following a path similar to the mode 10 locus 
as shown. A possible solution would be to move the actuator pole by redesign of the 
actuator system. This would require analysis of even higher frequency modes (assum- 
ing the actuator poles were pushed higher by redesign) to ensure that the trouble caused 
by actuator poles in close proximity to modal poles is not simply shifted to a new area 
of the s-plane. 

For the aft rate sensor location (Figure 52) all elastic mode loci exhibit stable 
configurations. Several modes do not have significant stability margins, but for this 
time of flight there are no problems. 

In order to complete the planar study, and to develop alternate sensor locations 
should the need arise when yaw and roll are realistically coupled, some combinations 
of sensors located on the orbiter and the tank were briefly investigated. One considered 
a combination of attitude sensor mounted forward on the orbiter (model node 22) and 
the rate sensor mounted forward of the first mode antinode on the tank (model node 18). 
The result was somewhat similar to the results for both sensors mounted forward on 
the orbiter (Figure 51). Here the orbiter first body mode (mode 10) again departs 
directly toward the right-half plane, but additionally modes 3, 5 and 8 also depart 
primarily to the right. The conclusion is that this sensor location combination has 
no advantages over the forward location on the orbiter. 

Another sensor location combination that was investigated had both rate and 
attitude sensors located forward on the tank at model node 18. Here modes 3 and 4 
depart to the right but have relatively small residues and could be gain stabilized. 
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Figure 51. Yaw Forward Rate Sensor (Node 14) and Forward Attitude Sensor (Node 22) on Orb iter 
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Figure 52. Yaw Aft Rate Sensor (Node 6) and Forward Attitude Sensor (Node 22) on Orbiter 





Modes 6 and 8, however, depart directly toward each other and have such large modal 
amplitudes that gain stabilizing them would be difficult. Any phase shift introduced to 
improve the stability of mode 6 will cause mode 5 to "swap" and depart into the right- 
half plane. Thus, this configuration also provides no apparent benefits. 

When the rate sensor was located aft on the tank (model node 5), and the attitude 

sensor mounted forward on the orbiter (model node 22), mode 10 presented the same 
problem as before; i. e. , the locus departs directly into the right-half plane with large 
modal gain. In addition the first, sixth, and seventh modes also depart somewhat to 
the right although they have smaller modal gains. 

Figure 53 shows the effects of locating both sensors aft on the orbiter (model node 
6). This configuration has results similar to those of Figure 52 . Although the margins 
are not large, at the required operating gains, all modes are stable without further 
compensation. 

Thus, using an aft rate sensor with either a forward or aft attitude sensor will 
provide a planar yaw (artificially decoupled) channel that is stable for the first twelve 
modes but does not have acceptable stability margins. 

Roll Planar Analysis . — The root locus diagrams for planar control of the roll 
channel are presented in Figure 54 for a forward rate sensor location (model node 22) 
and in Figure 55 for an aft rate sensor location (model node 6). The attitude sensor 
location for both of these cases is at model node 22 in the orbiter nose. 

In comparing these two figures, mode 7 shows the classic 180-degree phase 
shift associated with first body modes when the rate sensors are on opposite ends of 
the vehicle. Mode 9 shows a large modal amplitude with the rate sensor forward as 
compared to the small modal amplitude with the rate sensor aft. Mode seven in the 
aft rate sensor case is unstable at the closed-loop gain required for control. However, 
the aft sensor case appears more amenable to filtering because, in the forward rate 
sensor case, modes 7 and 9 depart directly toward each other, and the large modal 
amplitude associated with each will result in ’’loci swapping" as filtering is introduced 
to rotate the loci. With the aft sensor this does not occur, and Figure 56 shows the 
results of adding a quadratic lag filter to the aft rate sensor case. The filter has £= 0.7 
and co n = 117 rad/sec. This stabilizes the mode 7 locus to neutral stability but causes 
the locus associated with the actuator to rotate clockwise, causing the closed-loop 
operating point to lie in the right-half plane. Some further changes to the engine 
actuator system would be required to stabilize this configuration. 

Both attitude and rate sensors are located aft (model node 5 on the tank) as a 
possible alternate location (Figure 57). This is a satisfactory location, al- 
though modes 5 and 11 depart substantially to the right. The modal amplitudes for 
all modes are small, and at the gains required for closed-loop control all are 
slightly stable. 
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Figure 53. Yaw Rate and Attitude Sensors Aft on Orbiter (Node 6) 
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Figure 54 . Roll Rate and Attitude Sensors Rirward on Orbiter (Node 22) 
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Figure 55. Roll Rate Sensor Aft (Node 6) and Attitude Sensor Forward (Node 22) on Orbiter 
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Figure 56. Roll Rate Sensor Aft (Node 6) and Attitude Sensor Forward 
(Node 22) on Orbiter, Autopilot Compensation Included 



Another alternate combination that was investigated had the rate sensor aft on 
the tank (model node 5) and the attitude sensor aft on the orbiter (model node 6). This 
configuration had results similar to those obtained with both sensors aft on the tank. 
The mode 7 response is more pronounced but i departs to the left. Modes 5 and 11 
continue to depart to the right but again have operating points slightly in the left-half 
plane. 


A forward attitude sensor (model node 22) on the orbiter and an aft rate sensor 
(model node 5) on the tank has a pronounced mode 7 response departing directly into 
the right-half plane. Modes 5 and 11 have results similar to the other cases where the 
rate sensor was tank mounted. 

Finally, a case was considered in which both sensors were forward mounted, 
with the attitude sensor at model node 22 on the orbiter and the rate sensor at model 
node 13 on the tank. Modes 5 and 7 depart directly toward the right-half plane, and 
mode 7 is unstable at the required closed-loop gain. 

In the roll planar analysis, an aft rate sensor and an aft attitude sensor proved 
more satisfactory than the other conditions analyzed; however, none provided really 
satisfactory stability margins. 

Yaw/Roll Coupled Analysis . - The analysis for coupled yaw/ roll for this con- 
figuration follows substantially the same technique as that employed for the dual fly- 
back configuration. 

For an initial investigation, the roll and yaw forward attitude sensors were lo- 
cated at model node 22 on the orbiter and the aft rate sensors were located at model 
node 6 on the orbiter. This combination is considered as the most desirable from an 
equipment location point of view. Figure 58 shows the root locus diagram for this con- 
figuration with the roll gain (Kar) fixed at 10 and the yaw gain (Kay) varied. The effect 
is to move one mode loci completely into the right-half plane (at to = 26 rad/sec) and to 
have the closed-loop operating point of another (u) - 18 rad/sec) in the right-half plane. 
This is not a satisfactory sensor location autopilot gain combination. 

The second configuration investigated combined the aft yaw rate sensor (model 
node 6) with the forward yaw attitude sensor (model node 22) on the orbiter and aft 
roll rate (model node 5) with aft roll attitude (model node 5) on the tank. Figure 59 
shows the results of fixing =10 and varying Kay* 111 case the mode with 
poles at oo- 19 rad/sec has its operating point in the right-half plane, and the mode 
with poles at to = 28. 6 rad/sec has its operating point on the joo axis. 

In order to pinpoint just what was causing this coupled instability, Kar was set 
equal to zero and Kay varied. This removes the coupling through the roll autopilot. 
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The result is shown in Figure 60 and is very similar to that seen in Figure 59, but 

considerably different from the planar results shown in Figure 52, This would indicate 
that the problem is coupling through the yaw attitude control system. To ascertain this, 
Kay was set equal to zero and Kar varied. Figure 61 shows the results, which are 
similar to the corresponding planar analysis results (Figure 57). This verifies that 

the problem is yaw autopilot coupling. Since the yaw gain is primarily set at a high 
level due to the rigid-body restriction, there appears to be little chance of stabilizing 
this configuration with conventional autopilot techniques. 

No indications of whether a satisfactory system could be obtained with SRM TVC 
or with multiple sensor systems were obtained, since these configurations were not 
investigated. 

The marginal stability conditions illustrated in the planar analyses, the complexity 
of this coupled system, and the simplified nature of the elastic model indicated that 
further analysis at this coupled level was not warranted at this time. The results of 
this brief yaw/ roll coupled analysis indicated that modal coupling through the attitude 
control system is a very important consideration for this solid motor vehicle. 
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Figure 60. Coupled Yaw/Roll, Aft Yaw Rate Sensor on Orbiter, Roll Gain Set to Zero, 
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6 CONCLUSIONS 


Dual Flyback Vehicle 

The results of this study indicate that conventional autopilot techniques can be 
used to stabilize the vibration modes at the liftoff flight condition. Furthermore, use 
of the forward rate sensor configuration appears feasible in both the pitch and yaw 
autopilot control planes. This particular configuration has the greatest potential for 
later flight times when, if the elastic modes follow conventional patterns, the stability 
margins increase as modal frequencies increase. 

Propellant sloshing did not become an important stability consideration at this 
time of flight. No coupling is indicated between these modes and the lower frequency 
elastic modes. 

Elastic coupling between the orbiter and the booster does not appear to appre- 
ciably affect the total vehicle elastic mode stability. ' This is probably due in large 
part to the fact that the orbiter weight is only 20 percent of the booster weight at lift- 
off. Thus at later flight times, when the vehicle weights are more comparable, signif- 
icant changes in the elastic coupling between the two vehicles could occur. This obser- 
vation and the omission of aeroelastic effects means that the results of this study can- 
not be extrapolated with any degree of confidence beyond the liftoff flight condition. 


Solid Motor Vehicle 


Due to the preliminary nature of the modal data available and the extreme sensi- 
tivity of control system stability to these data, the results of this study are inconclu- 
sive in determining whether or not conventional autopilot techniques can be used to 
adequately stabilize the vehicle. 

Although positive elastic mode stability is attained for the symmetric mode case, 
the margins of stability are small. For the antisymmetric mode case significant 
elastic mode coupling is evidenced for the yaw/roll coupled system. This coupling pro- 
duces unstable modes, which are not amenable to conventional stabilization techniques. 
Propellant sloshing is not an important stability factor at liftoff. 

Three important interrelated considerations that complicate the elastic mode 
stability of this vehicle are the high modal gain, the location of the "tail-wags-dog" 
zero, and the high degree of elastic coupling between the vehicle elements through the 
links that join these elements together. The high modal gain is due in large part to 
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the relatively high autopilot forward-loop gain required to obtain satisfactory control 

mode performance characteristics. This high autopilot gain is, in turn, due to the 
low control thrust (considering the size of the vehicle) provided by the orbiter engines. 

The "tail-wags-dog" zero is located in the frequency range of the lower fre- 
quency body modes. The phase shift introduced by this singularity not only compli- 
cates the stability problem for the liftoff condition, but it will also introduce extra- 
ordinary considerations as the modal frequencies change with increasing flight times. 

The combination of the high modal gain and the "tail-wags-dog" zero location 
accentuates the importance of the elastic mode frequency density, increases the 
potential coupling through the attitude control system, and increases the sensitivity 
of the system stability problem to variations in the modal parameters. 

Based on these approximate mode shapes, and the preceding comments 
conventional, single-rate-sensor stabilization does not appear feasible if only the 
orbiter engines are used for control. However, a brief investigation of adding solid 
motor TVC (secondary injection) indicates that improved stability properties can be 
attained with this increase in control thrust. Another brief study of a two rate sensor 
autopilot configuration indicates that the vehicle stability properties can be improved 
by employing a multiple sensor autopilot configuration. 
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7 RECOMMENDATIONS 


The complexity of the stability problem for the solid motor vehicle and the 
choice of this concept for continued development indicate that studies of vehicle vi- 
bration mode stability should be continued. These studies should be expanded to 
include other critical flight times (e.g. maximum dynamic pressure, solid motor 
burnout, etc . ) . They should also include alternate autopilot configurations, such as 
multiple sensors and adaptive techniques. 

The brief investigation of multiple sensors showed sufficient improvement that 
further study of this concept is definitely warranted. The use of multiple sensors 
for redundancy has been considered in earlier space shuttle studies, although all 
sensors were at the same location. This concept should now be broadened to a study 
of multiple sensors at various locations to satisfy both redundancy and elastic mode 
stability requirements. 

The structural complexity of this solid motor vehicle concept and the sensitivity 
of the vehicle stability to its elastic properties indicate that further elastic model 
development should be pursued in parallel with the stability studies. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Virginia, November 24, 1972 
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APPENDIX A 


Symbols 
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LIST OF SYMBOLS 

Definition 

Gravity constant 

Engine gimballed moment of inertia about gimbal point 

Vehicle roll moment of inertia about center of mass 

Vehicle pitch moment of inertia about center of mass 

Vehicle yaw moment of inertia about center of mass 

Vehicle x-z plane product of inertia about center of mass 

Autopilot forward loop gain 

Actuator servo loop gain 

Actuator load torque feedback gain 

Autopilot rate -to -attitude gain ratio 

Attitude sensor gain 

Rate sensor gain 

Autopilot forward loop static gain 

Sloshing pendulum length 

Aerodynamic roll moment coefficient due to perturbation of 

Distance from vehicle center of mass to liquid control engine 
gimbal point along X body axis 

Distance from vehicle center of mass to solid engine thrust 
application point along X body axis 
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APPENDIX A — Continued 


l Ez 

^Ey 

SL-p 

% 

^RE 

m e 

M p 

Mr 

m t 

M u 

% 

% 

N /3 

s 

T (or T c ) 
T RAT 


Distance from vehicle centerline to liquid control engine gimbal point 
along Z body axis 

Distance from vehicle centerline to solid engine thrust application 
point along Y body axis 

Sloshing pendulum moment arm 

Distance from gimbal point to engine center of mass 

Equivalent roll control moment arm 

Engine gimballed mass 

Propellant sloshing mass 

Vehicle reduced mass (M^ - 2M p ) 

Vehicle total mass 

Aerodynamic pitch moment coefficient due to perturbation of u 
Aerodynamic pitch moment coefficient due to perturbation of w 
Generalized mass of ith elastic mode 
Vehicle steady-state pitch rate 

Aerodynamic yaw moment coefficient due to perturbation of (3 
Laplace variable (complex frequency) 

Liquid engine control thrust 

Ratio of equivalent roll control thrust to total liquid engine control 
thrust 

Solid engine thrust 

Vehicle steady-state forward velocity 
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APPENDIX A — Continued 


l L 


Y rp 




X, 


u 


X, 


w 


Z u 

Z w 

Z cg 

6 

op 

5 oy 

£cn 

6 n 


^PA, %A 


a RA 

a PRS, 0r yRS 
O' 

RRS 


Vehicle steady-state forward acceleration 
Magnitude of vehicle velocity vector 
Vehicle steady-state plunging velocity 

Aerodynamic longitudinal force coefficient due to perturbation of u 

Aerodynamic longitudinal force coefficient due to perturbation of w 

Aerodynamic side force coefficient due to perturbation of £ 

Aerodynamic plunging force coefficient due to perturbation of u 

Aerodynamic plunging force coefficient due to perturbation of w 

Distance from vehicle centerline to vehicle center of mass along Z axis 

Steady-state thrust vector angular deflection in pitch 

Steady-state thrust vector angular deflection in yaw 

Actuator servo damping ratio 

Actuator system closed loop damping ratio 

Steady-state pitch angle relative to launch pad 

Density of liquid propellants 

Elastic mode slopes at attitude sensor locations in pitch and yaw 
control planes 

Elastic mode roll deflection at roll attitude sensor location 

Elastic mode slopes at rate sensor locations in pitch and yaw 
control planes 

Elastic mode roll deflection at roll rate sensor location 
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APPENDIX A ~ Concluded 

Elastic mode slopes at thrust application point in XZ and XY planes 

Elastic mode roll deflection at thrust application point 

Elastic mode deflection at propellant sloshing hinge point parallel 
to vehicle Z and Y ax es 

Elastic mode deflection at thrust application point parallel to vehicle 
Z , Y and X axes 

Elastic and propellant sloshing natural frequencies 
Actuator servo resonant frequency 
Frequency of tail-wags-dog zero 
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APPENDIX B 

ANALYSIS MODELS AND PERTURBATION EQUATIONS 

This appendix summarizes the analysis models and presents the perturbation 
equations that are used in this analysis. It also briefly describes the solution tech- 
nique that is used to calculate the roots of the characteristic equation. 

The definition of the vehicle coordinate system is shown in Figure Bl. The 
origin of the coordinate axes is located at the vehicle center of mass. 

RIGID BODY EQUATIONS 
Translation Equations 

AF 

su+Q w+W s 0 = - X - 
o o 

T 



z 


k b 

Figure Bl. Vehicle Coordinate System 
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sw - Q u - U s0 
o o 


AF 

z 


M 


R 


A F = S (thrust forces) + E (gravity forces) 

+ E (aerodynamic forces) + E (engine 
inertia forces) + E (slosh forces) 

Rotation Equations 


2 xz , 2 , z z vv 

S 0 - y~ (s 0 + Q q S0) + — — Q S ip 

XX XX 


A L 


XX 


„ AL 

2 V 

s e =-y JL 

yy 


9 I o I - I AL 

.-22 (s 2 * - q s4) > + J2L_JH5 Q s<t ,=_^ 

1 0 I o I 

ZZ zz zz 


AL =E (thrust torques) +E (aerodynamic 
torques) + E (engine inertia torques) 

+ E (slosh torques) 

The engine inertia forces and torques used in these equations are only those 
due to engine gimbal angles. 


In the above equations: 



U = U 0 + u 

P = Po + P 

p =0 

W = W 0 + w 

Q = Qo + Q 

q = e 

ii 

< 

c 

+ 

< 

R = R 0 + r 

r = $) 
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Longitudinal Equation 


2 2 1 
(s + 2 £ CO S + 60 ) F , = -- 8 W 

pk pk pk pk L 

pk 


f & - L U 

pk pk 2 o 


+ ' s + ~ s'l 6 -S -1 <i>. s 2 q. 

L , L , / j lzpk i 


pk 


pk 


Lateral Equation 


,2 , 2 U o 

,s + 2 f P k“pk B + a, pk )r pk = 'ir 


/H - L U \ 

r pk pk 2 o', 

s + — 

V pk pk • i "pk 


pk 

'12 
s ) $ - £ rr— s q. 

/ i , iypk \ 


2 U o 

In these equations to , = — - and k = 1, 2. . . . , m 

pk L 

pk 


ELASTIC VEHICLE EQUATIONS 


2 2 

(s + 2 £ W s + 0) ) q = — 

^ Bi Bi Bi ' M i Tf]_. 


Q. = £(thrust forces) • 


+ £ (engine inertia forces) • <i>. 


Ei 


+ £ (sloshing forces) • $ 


£ (engine inertia torques) * a 
3<i>. 


Ti 


where a. - —rr and i - 1, 2, . „ . , n 

1 a A 
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APPENDIX B — Continued 


The engine inertia forces and torques used in these equations are only those due 
to engine gimbal angles. The modal displacements and slopes, cp.-ando - ., are obtained 
from orthogonal normal modes. 

GRAVITY PERTURBATION FORCES 



Z LAUNCH SITE 

COORDINATE SYSTEM 


F = - M T g (cos 0 ) 0 
gx 1 o 

F = M T g [(sin 0 )0 + (cos 0 ) 0] 
gy 1 o o 


F = - g (sin 0 ) 0 
gz T o 


Since the center of the body axis coordinate system is located at the vehicle 
center-of-mass the gravity torques equal zero. 
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APPENDIX B - Continued 


AERODYNAMIC PERTURBATION ACCELERATION COEFFICIENTS 


The aerodynamic forces and moments are expressed as dimensionalized coeffi- 
cients (l/sec) for the body variable rates (u, w, 9, etc.). 

Translation Equations 


F 

— = X 
M t u 


u + X w 
w 


77T7 = (Y • s + Y ) /3 + Y s0 + Y s >p 
U q M r p j9 M P r 


F 

za 

- — = Z u+(Z.s + Z)w + Zs0 
M r u w w 0 


Rotation Equations 


xa 


~ — - Lfl+L S0 + L s 4> 
I R H pa Y r 

xx ^ 


~~~ = M u + (M. s+M ) w + M s 0 
I u w w 0 

yy 



= (N. s + N)S+N S0 + N sill 
0 0 P P r v 


ENGINE INERTIA PERTURBATION FORCES AND TORQUES 

Figures B3, B4, and B5 illustrate the engine gimballing dynamic models and 
parameter definitions. 
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(a) Pitch and Yaw Control 




<L = UNDISTORTED BOOSTER CENTERLINE 

<|/ = DISTORTED BOOSTER CENTERLINE ,, 

4 = £ $(x, z)jrp qj — 

1 k b 

(Z) 


Figure B4„ Longitudinal (Pitch Plane) Rocket Engine Dynamics 





APPENDIX B - Continued 


S$ixT Qi 


i b < L 
(X) 



(J, = UNDISTORTED BOOSTER CENTERLINE 
qf = DISTORTED BOOSTER CENTERLINE 

l = L 4>(x, y) iT qi 
i 


Figure B5. Lateral (Yaw/Roll Plane) Rocket Engine Dynamics 
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The engine inertia forces and torques acting on the vehicle at the gimbal 
point are as follows: 

F = - M su - [m (Z - z )-M J i h ] s 2 fl 
Ex E E ' Ez eg' E R°op 


S[M Z 6 a. + M 1 s 2 q 

j E R op izT E ixT^ i 


-M Z b s 6 
E R op p 


F — M U + M (& - z - z 6 ) s 0 

Ey E o E v Ez eg R op' v 


+ C~ M (Z ~ Z ) s -MU s+M U]ii) 
E C R T7 ^ TP T 


E o 


E o J 


+ £ [ -M <£. + M z a. _ - M_ Z _ 6 a ] s q 

i E lyT E R lyT E R op ixT 4 i 


M_ Z s 6 
E R y 


F Ez = - M E sw+[M E<VV B +m e u o s 


-M U ] 0 +£ [ - M $. + M Z a . J s q 

E o i '- E izT E R izT M i 


+ M Z s 6 
E R p 


, =-M U Z 6 sfl+C-L + M n jL6 (A„-z )] s 0 

Ex E o R op r R ER op Ez eg 

+ [-M Z z 6 s 2 - M Z n U 6 s + M_ Z u 6 ] ih 
E R c op E R o op ' E R o op y 


+ £ [-M Z 6 $ - I g ] s q 

i E R op iyT R u ixT J 4 i 


+ T jyr Z Z s 6 

RAT E R RE r 


116 



APPENDIX B - Continued 


L = M & 6 s u - M & s w + [ M L 6 (£ -z ) 

Ey E R op E R E R op Ez eg 

+ - I ]s 2 0+B{[ - M 


E R c R 


E R ~izT 
2 


+ L CT . M fi A. J S 

R izT E R u op nxT 

+ M E ^ u 0 °i Z T s + M E i a izT } «i 

+ (lRS 2 -M EiR U oS+ M E . R 6 o ,6 p 

L E Z =M E U o^ S ? + [- M E\^ z -V s2 
+ M E i R 6 cp U o S - M E‘ R «„p U o ]0 
""VrW =% + S[[M E £ B $. yT 
2 

-I_, a. ] s + [-M $. U a. + M j0 6 U ct ] s 

R iyT E R o iyT E R op o ixT 

+ [ - m e *r *0 V- + m e *r °o 6 op °ixT ] 1 q i 

+ ri K s 2 + M E ^ u oS + M E i K u o ]6 y 

The total engine inertia torques acting on the vehicle are obtained from the 
following expressions: 

L = L - (£ - z ) F 

Rx Ex Ez eg Ey 

L_ = L + (i ~z ) F - H F - E M U <&. q. 

Ry Ey Ez eg Ex c Ez 1 E o xzT l 


L Rz L Ez + l c F Ey 
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THRUST PERTURBATION FORCES AND TORQUES 
The thrust forces and torques are obtained from Figures B3 S B4, and B5„ 

F =-LT6 a. q. - T6 6 
Tx • op xzT 1 op p 


=£(Ta. m - T 6 a. m + T cr. m ) q 
Ty . lyT op ixT s lyTs i 


- T6 


F = S (T a. _ + T n m ) q. + T 6 
Tz . izT s izTs l p 


L [Tfi $ - T (£ -z ) a 

Tx j op lyT Ez cg ; lyT 


T (A - z ) 6 a. ^ + T 6 

kz eg op ixT s oy l 


oy iz Ts 


+ T L a. + T z a. m 
s Ey izTs s eg lyTs 


"T i 6 a. ] q. + T (i - z ) 6 
s Ey oy ixTs l Ez eg y 


+ T RAT T £ RE § r 


L_, = E[T$. -T6 (jL -z ) a. m 

Ty . izT op Ez eg izT 


-Tj i a. " T 6 § + T 

c izT op ixT s izTs 


™T i a. _ ] q. + [T 6 (z - l ) 
s c izTs l op eg Ez 


-T i ]6 
c p 
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Tz 


£[-T4 +TU. "T| fi CT. 


lyT 


c iyT 


e op ixT 


T <£> + T (j l - JL„ 6 )CT. m 

s lyTs sc Ey oy lyTsj q. 


- T t 6 

c y 


PROPELLANT SLOSHING PERTURBATION FORCES AND TORQUES 
The propellant sloshing forces and torques are obtained from Figure B2. 


F =0 
Sx 


F - U £ M T , 
Sy o k pk pk 


F = U £ M r 

Sz o k pk pk 


l = u s m z r . 

Sx o k pk eg pE 


L =-U £ M L T. 
Sy o k pk pk pk 


L = U £ M IT. 

Sz o k pk pk pk 

ACTUATOR PERTURBATION EQUATIONS 


The same third-order linear actuator model is used in all three control planes. 
Figure B6 illustrates this model and presents the parameter values used in this study. 
The log magnitude and phase diagrams for this particular set of parameters are given 
in Figure B7. 
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APPENDIX B - Continued 



co ^ = 50Rad./Sec. (Closed-Loop Resonant Frequency, at Load) 

C d = 0.3 (By Design) 

A =25 (1/2 u! ) 
d 

o) = 70 Rad/Sec (Servo Resonant Frequency) 

= 12. 8 Sec -1 - (Servo Loop Gain) 

K q = 55 Sec * (2 C d w d + A) 

Figure B6. Actuator Dynamic Model 
The perturbation equation is : 


3 2 2 2 
(s + 2 C to s + oo s + K oo ) 6 
cn cn cn c c 


= K oo 2 6 + (s + K )-f- T (s) 
c c c oIL 

R 

where* 2 C oo = 2 C . £0 + A 
cn cn d d 

2 2 

oo -to + 2 £ co A 
cn d d d 


9 9 

K oo = A oo 2 
c c d 


An alternate actuator model uses only a single-order lag: (s + K ) 6 = K c 6 C - 
This model is adequate for low frequency modes, such as rigid body and, sometimes, 
propellant sloshing. 
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APPENDIX B — Continued 


The load torque feedback for each control plane is given below. 

Roll Control Plane 


T LR ? M E £ R ^RE °ixT S Q i 


Pitch Control Plane 


T„ ^ - M„ ^ 6 s u + M„ s w 


LP 


E ~R op 


E 




+ [L - M_ Jt - M„ e 6 tf _ - 2 )] s 6 

R E R c E R op Ez eg 


(- 1 

+ X /M l 3> - I ct -M 1 6 Js q. 

f U E ^R “ zT R izT E R op ixT 


Yaw Control Plane 


T LY = - M E *R U o s ^ M E 'b “ez ■ V s 4 


- ,m e fa ' V s * 


+ 


f <- M E ^ V + b "lyT* S q i 


CONVENTIONAL AUTOPILOT 

The equations that describe the conventional autopilot are based on the block 
diagram shown in Figure B8. Although the following equations were written specifically 
for the pitch control channel, these equations are also applicable to the other two con- 
trol channels when the appropriate changes are made to the perturbation variables. 

Rate Sensor 


1 (*1 

ir Cr i <s)]e Rs =s0_ f CT ms sq i 

R l 

(*) For the yaw and roll control channels this sign is + based on the sign con- 
ventions shown in Figures B4 and B5. 
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APPENDIX B — Continued 


Attitude Sensor 
[F 2 (s)]e A =9 < : ) La tA c, i 


(*) For the yaw and roll control channels this sign is + based on the sign con- 
ventions shown in Figures B4 and B5. 

Actuator Command 


CF 3 <S)]6 c = - K A e RS- K A 0 A 



Figure B8„ Attitude Control System. Block Diagram 
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APPENDIX B - Continued 

In the preceding equations, F 1 (s), F 2 (s) and F 3 (s) indicate capability for the 
addition of compensation filters. 

SOLID MOTOR THRUST CONTROL 

Secondary injection is considered for solid motor thrust control in the symmetric 
mode case (longitudinal equations) . To incorporate this control function into the sys- 
tem equations requires the following additions to the longitudinal equations. 

w Equation 

F = T 6 
z 6 TS S PS 

0 Equation 

L y 6 TS = ' T S *cS 6 PS 
q, Equation 

Q = T $ 6 

16 TS S izTS PS 

The TVC system is approximated by: 

< s + V *P S = K o 6 CP 

where K = 15 sec (t = 0.0067 sec) 
c 

The command for the TVC system is the same as that sent to the liquid pro- 
pellant gimballed engines. 

COMPLETE SYSTEM EQUATIONS 

The complete system equations, as they are used in this study, are presented 
here. Terms which were demonstrated to have a negligible effect on the analysis re- 
sults have been eliminated. The primary reason for this simplification is to reduce 
the complexity and cost of obtaining the roots to the characteristic equation. 
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These equations can be represented in a matrix format. The roots of the 
characteristic equation are determined by finding the eigenvalues of the determinant 
on the left-hand side. The root loci are then obtained by plotting these roots, for 
different values of autopilot gain, on the complex plane and connecting the points in 
the direction of increasing gain values from 0 (open-loop pole) to » (zero). 


LONGITUDINAL EQUATIONS 


u Equation 


M e 6 od 2 

(s - X ) u + (Q - X ) w +\ r r~ s + W s + g cos 0 10 

u o w \ o o/ 


+ £ 
i 


f M 4 6 \ 9 /T6 

■5-*L°Eor ■ s 2 + __2I> 

M„ 12 1 / \M„ / izT 


M l 6 0 T 6 

, , E R op 2 op , . 

q +i s + ^ 6 =0 

i V M t M t / °p 


w Equation 


(- Q + Z ) u + (s - Z ) w + (- U s + g sin 0 ) 0 
o u w o o 


° / T T 

M ^ M pk^pk f ( M °izT M °izTS ) ^i 
H * \ R R / 


M F^R 2 T 

+ 1 5_5. s 2 -.— 16 = 0 


M. 


R 


m r /P 


0 Equation 


M„ 6 


EX op 

— 7 s-M Ju-M w + 

I u / w 

yy 


1-2 


M A_6 (l_ - Z > 

E R op \ Ez eg; 

I 

yy 


2 

s 0 


U 

+ ~T~ EM. i , r . + E 

J yyk Pk Pk Pk i 


MJJ (A -z ) M A 6 

E R op Ez eg E R op, 

— — rr $ 

I izT I ix T 


yy 


yy 
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T l T 6 (A_ - z )\ T 6 

4 +1—5+ — 2E_^_££ L + — 2E * 

I izT \I I / izT I ix T 

yy \ yy yy / yy 


T S , + T S *CS 

OP + cr 

I lzTS I is 

yy yy 






q. + 



i 



m i i i m a 6 (A_ - z ) 

E c R _R E R op Ez eg 

I 1 I 

yy yy yy 


M -^R U o 


E 


s + 


yy 


t i T 6 - z ] m„ u 


c op\ Ez eg j E R o 


yy 


yy 


yy 


6 - 0 
P 


Fpk Equation (Propellant Sloshing) 


Ju_ s w + ( - r — s 2 - ~ s )e + (s 2 + 2^ co s + co 2 ) F 

V V V V / pk pk pk pk 


1 2 

+ S r — <£>. s q. - 0 

. L , lzpk i 
i pk 

2 * 0 

where co , = - — and k = 1, 2, . . . , m 

pk L 

pk 


q. Equation (Elastic Body) 


M 


„ iL 6 cr. s u + [M_ H_ 6 (i - z ) a. - M_ i, 6 4>. ] s 0 

E Tl op izT L E op v Ez eg' izT E 1R op ixT J 


- 6 0 S k M pk r pk + % (B + 2 t Bi “Bi S + ' “ Bi> <4 

n 2 

C (M E 4. V *ixT> %T + (M E 4 s 6 op°izT> V’ S 

+ [(- T «S>. zT + T 6 Qp ^.j.) ct j zT “ ( T s ®i zT g) °jzTS^ q j 


+ [(“ M_ £ _ + I_ a. + M_ JL 6 _) s + M_, L U cr. s 

E R izT R izT E 11 op ixT E K o izT 
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+ (- T $ + M i U a + T 6 J] 6 = 0 

izT E R o izT op ixT p 


5 Equation 

P 


M_ JL 6 _ M $, 

E R op , 2 , E R 2 

M E Vc 

(s + K s) u - ■■■ ~(s + K s) w + 

L 0 I 0 

R R 

— ~ 1 
'r 


M a 6 (jL, ~ z ) 
E R op Ez eg 


I. 


R 


3 2 

(s + K s ) 0 + 2 

° i 


M„ A 

E R _ 

— - <P. m + a. 

I izT izT 
H 


M A 6 

+ ■ E R op $ 

I ixT 

XX 


3 2 3 2 2 2 

(s + K s ) q. +(s +2£ CO S+co s + K co ) 6 
o l cn cn cn c c p 


K co 6 =0 

c c cp 


PITCH AUTOPILOT EQUATIONS 
Rate Sensor 


1 / 1 2 

1F~ ( — ~9 s + s + 1 j 0 - s 0 + E <7 s q. = 0 

K \ 2 60 Til / RS • lPRS 1 

RP\co FI / 1 


Attitude Sensor 
1 2 2 ^F2 \ 

~ s + s + 1 ]0 - 0 +£ CT q. = 0 

2 C0j_„ / A . iPA A 

W F2 / 


Actuator Command 


1 2 2 ^F3 

g + - - s + 1 ) 6 + K 0 +K 0=0 

2 m C p p RS p A 

„ 60 f3 f3 
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LATERAL EQUATIONS 


|3 Equation 
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0 Equation 
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y 


0 Equation 
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T pk Equation (Propellant Sloshing) 
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q. Equation (Elastic Body) 
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ROLL AUTOPILOT EQUATIONS 
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APPENDIX C 


DUAL FLYBACK VEHICLE DATA 


This appendix contains a tabulation of the dual flyback vehicle data used in this 
study. It also contains illustrations of the pi’opeilant sloshing tank geometries, the 
specific elastic vehicle models, and the elastic mode shapes. 

Essentially all of these data were generated during the Space Shuttle Phase B and 
follow-on studies. The customary (English) system of units was used to develop the 
data. 

The definition of symbols is contained in Appendix A. . 

The propellant sloshing parameters were generated by a version of the digital 
computer program described in Reference 2. 

The symmetric configuration uses a flexible beam model, whereas the anti- 
symmetric configuration uses a detailed three dimensional model for the booster 
vehicle. Both configurations use a simplified eleven-node flexible beam model for 
the orbiter vehicle. The mode shapes were generated by the NASTRAN computer 
program. Further discussion of these models and of the NASTRAN computer program 
is contained in Reference 3. 

TABLE Cl . - VEHICLE MASS PROPERTIES 


Parameter 

Value (Liftoff + 10 seconds) 

T (12 engines) 

2.94 X 10? N 

(6.61 X 106 lb) 

Mg (12 engines) 

4.48 x 10 4 kg 

(3070 slugs) 

m R 

2 . 06 X 10 6 kg 

(1.412 x 10 5 slugs) 

M'p 

2.21 X 10 6 kg 

(1.515 xlO 5 slugs) 

I-^ (12 engines) 

5. 86 X 10 4 kg-m 2 

(4.32 X 10 4 slug - ft 2 ) 

■*-xx 

4.53 x 10 7 kg-m 2 

(3.34 X 107 slug - ft 2 ) 

St 

7.32 x 10 8 kg-m 2 

(5.40 x 10 8 slug - ft 2 ) 

*zz 

7 . 16 X 10 8 kg-m 2 

(5.28 X 10 8 slug - ft 2 ) 

Scz 

-1.970 x 10 7 kg-m 2 

(-1.453 X 10? slug - ft 2 ) 

*c 

-44,5 m 

(-146 ft) 

A Ez 

0 m 

(Oft) 

£ r 

0. 826 m 

(2.71 ft) 

%E 

4. 61 m 

(15.14 ft) 

z cg 

-1,22 m 

(-4 ft) 

6 op 

-0. 07 rad 

(-4 deg) 

t rat 

2 
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APPENDIX C - Continued 


TABLE C2 . - STEADY-STATE TRAJECTORY DATA 


Parameter 

Value (Liftoff +10 Seconds) 

g 

9.81 m/sec 2 

(32.17 ft/sec 2 ) 

Qo 

-0.0035 rad/sec 

(-0.2 deg/sec) 

Vo 

31.09 m/sec 

(102 ft/sec) 

U 0 

13.14 m/sec 2 

(43.1 ft/sec 2 ) 

v T 

31.09 m/sec 

( 102 ft/sec) 

W o 

1 . 86 m/sec 

(6.1 ft/sec) 

e o 

1 . 51 rad 

(86.8 deg) 

MACH 

0.09 


Altitude 

151 . 8 m 

(498 ft) 

Dynamic Press. 

574.6 N/m 2 

(12 lb/ft 2 ) 

TABLE C3 , 

- TOTAL VEHICLE AERODYNAMIC DATA 

Parameter 

Value (Liftoff +10 seconds) 

X u 

-0. 0006 1/sec 


Xw 

0. 00012 l/sec 


Z u 

0. 0041 l/sec 


Z w 

0.144 l/sec 


Mu 

-0.000256 — — 

(-0.000078 — — ) 


m-sec 

\ ft-sec' 

M w 

0.00699 1 

( 0.00213 - 1 ) 


m-sec 

V ft-sec/ 

Y /5 

-0.430 m/sec 2 

(-1.41 ft/sec 2 ) 

L/3 

-0. 057 l/sec 2 


Nj3 

0. 044 l/sec 2 
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1339 S. N. 


1479 S.N. 


9. 83 M 
(32. 25 FT) 


1866 S.N. 


2006 S.N. 


NOTES: 

(1) ALL BULKHEADS ARE ELLIPTICAL SURFACES. MAJOR RADIUS =/2~ X MINOR RADIUS. 

(2) QUIESCENT FLUID LEVEL- USED IN STUDY (-> ~). 


2042 S.N. 


2182 S.N. 


38. 08 M 
(124.92 FT) 


3681 S.N. 


3821 S.N. 


Figure Cl . Main Propellant Tank Geometry 
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APPENDIX C - Continued 


TABLE C4 , - PROPELLANT SLOSHING DATA 


Propellant Tank 

Parameter 


Value 

Oxidizer 

M p 

1.51 X 10 5 kg 

(10 350.0 slugs) 

(L° 2 ) 

ir 

Lp 

1.98 m 

(6.5 ft) 



15.67 m 

(51.4 ft) 


Op 

2.58 rad/sec 

(0.411 Hz) 


P 

1136 kg/m 3 

(70.9 lb/ft 3 ) 

Fuel 

Mp 

12 565. kg 

(861.0 slugs) 

(LH 2 ) 

Lp 

2.71 m 

(8.9 ft) 



-4.30 m 

(-14.1 ft) 


COp 

2.2 rad/sec 

(0.350 Hz) 


P 

70,16 kg/m 3 

(4.38 lb /ft 3 ) 


TABLE C5 . - AUTOPILOT AND ACTUATOR DATA 

Parameter 

Value 

K c 

12.8 1/sec 

Ko 

55 1/sec 

U>c 

70 rad/sec 

w cn 

57 rad/sec 

^cn 

0.48 (N.D.) 

% 

0.7 sec 

k a 

variable 
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APPENDIX C — Continued 


LENGTH (m) 


]0 


20 


30 


40 


50 


60 



1,000 


2,000 


BOOSTER 

STATION 


3, 000 


4, 000 L- 


ORBITER 



Figure C2. Elastic Model Schematic (Symmetric Boundary Conditions) 









APPENDIX C - Continued 


TABLE C6 . - ELASTIC VEHICLE MODAL DATA 

(SYMMETRIC BOUNDARY CONDITIONS) 


Mode 

Number 

rri 

kg 

(slugs) 

CO 

rad /sec 
(Hz) 

$xT 
m . 

m* 10 

$ zT 

HI , „9 

- x i(r 

m 

$ zP 
m „ o 
— X 10 3 
m 

ct zT 

rad/m x 10 4 
(rad /ft x 10 4 ) 


1 

350.3 

7.2 

1.893 

-1.055 

2.841 

3.97 


(24) 

(1.15) 




(1.21) 

2 

350.3 

11.6 

-1.939 

2.631 

11.70 

-14.76 


(24) 

(1.85) 




(-4.50) 

3 

350.3 

14.6 

-2.769 

-1.234 

-9.128 

8.24 


(24) 

(2.32) 




(2.51) 

4 

350.3 

17.9 

-0.1641 

3.046 

-0.798 

-14.21 


(24) 

(2 . 85) 




(-4.33) 

5 

350.3 

20.1 

-0.029 

-0.7294 

1.303 

7.94 


(24) 

(3.20) 




(2.42) 

6 

350.3 

24.4 

-16.87 

0.1503 

-1.59 

-1.76 


(24) 

(3.88) 




(-0.536) 

7 

350.3 

26.6 

2.635 

0.9535 

-7.501 

-7.71 


(24) 

(4.24) 




(-2.35) 

8 

350.3 

32.4 

0.426 

0.2033 

0.977 

-2.87 


(24) 

(5.16) 




(-0.874) 

9 

350.3 

37.9 

1.655 

-3.03 

4.309 

37.96 


(24) 

(6.03) 




(11.57) 
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APPENDIX C - Continued 


TABLE C7. - PITCH CONTROL PLANE AUTOPILOT SENSOR ELASTIC 

MODAL DATA (SYMMETRIC BOUNDARY CONDITIONS) 


Mode 


°PRS or °pa 


rad „ a /rad 

it * 10 (ir 


Model Node 201 Model Node 204 


X 10 4 j 

Model Node 207 


Number S.N. 1000 S.N. 1860 S.N. 3000 


1 

0.745 

0.341 

2.920 


(0.227) 

(0.104) 

(0.890) 

2 

17.06 

12.47 

-8.924 


(5.20) 

(3.80) 

(-2.72) 

3 

-15.36 

-9.646 

4.593 


(-4.68) 

(-2.94) 

(1.40) 

4 

-1.125 

-0.3806 

-4.068 


(-0.343) 

(-0.116) 

(-1.24) 

5 

-3.281 

-4.528 

3.773 


(-1.00) 

(-1.38) 

(1.15) 

6 

-4.26 

-1.345 

0.219 


(-1.30) 

(-0.410) 

(0.0667) 

7 

-24.84 

-7.973 

0.951 


(-7.57) 

(-2.43) 

(0.290) 

8 

7.678 

3.806 

-1.207 


(2 . 34) 

(1.16) 

(-0.368) 

9 

20.11 

-0.5906 

-1.955 


(6.13) 

(-0.180) 

(-0.596) 
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MODE 4 THRUST STRUCTURE MODE oj = J7. 9 RAD/SEC (2. 85 Hz) 




APPENDIX C — Continued 
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MODE 7 ORBITERMODE oj = 26. 6 RAD/SEC (4. 24 Hz) 

Figure C3. — Continued 




MODE 9 SECOND BOOSTER BODY MODE gj = 37.9 RAD/SEC (6. 03 Hz) 



(a) Three Dimension Booster With Orbiter Flexible Beam Model 


LENGTH STATION 

im NUMBER i 



(b) Booster Centerline Representation (c) Fuselage Model Cross Section 

Figure C4„ Elastic Model Schematics (Antisymmetric Boundary Conditions) 
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APPENDIX C - Continued 


TABLE C8. - ELASTIC VEHICLE MODAL DATA (ANTISYMMETRIC 
BOUNDARY CONDITIONS) 


Mode 

Number 

rri 

kg 

(slugs) 

CO 

rad/sec 

(Hz) 

$yT 
m q 

— X 103 
m 

m q 

— X lO 01 
m 

°xT 

rad/m x 10 4 
(rad/ft x 10 4 ) 

a yT . 
rad/m x 10 4 

(rad/ft X 10 4 ) 

1 

350.3 

5.4 

0.4114 

1.63 

-21.52 

-0.275 


(24) 

(0.86) 



(-6.56) 

(-0.0837) 

2 

350.3 

7.4 

-11.59 

-6.79 

10.70 

4.03 


(24) 

(1.2) 



(3.26) 

(1.228) 

3 

350.3 

10.8 

22.25 

2.82 

3.81 

-13.35 


(24) 

(1.7) 



(1.162) 

(-4.07) 

4 

350.3 

11.7 

-8.327 

-8.99 

-17.16 

6.39 


(24) 

(1.9) 



(-5.23) 

(1.948) 

5 

350.3 

14.8 

-8.692 

-3.22 

-19.19 

3.895 


(24) 

(2.4) 



(-5; 85) 

(1.187) 

6 

350.3 

18.2 

16.72 

-4.29 

6.12 

-14.8 


(24) 

(2.9) 



(1.867) 

(-4.51) 

7 

350.3 

21.6 

-4.084 

-3.06 

23.39 

4.666 


(24) 

(3.4) 



(7.13) 

(1.422) 

8 

350.3 

22.5 

9.437 

2.04 

-5.73 

-9.269 


(24) 

(3.6) 



(-1.747) 

(-2.825) 

9 

350.3 

25.6 

-6.216 

0.935 

-72.5 

7.185 


(24) 

(4.1) 



(-22.1) 

(2.19) 


145 



APPENDIX C - Concluded 


x 


b 

rt 

U 


rH 

o 


o 

o 

rH 

0 

CO 

rH 

05 

CO 


© 


cq 

<N 

O 




CD 

tH 



CH 

t? 

O 


b 

C\| 

© 

05 

CD 

f=H 

CD 

cO 

co 

05 

LO 

oq 

CD 

05 

co 

CO 

rH 

rH 


o 

CO 

LO 

LO 

CD 

CM 

00 

rH 

r=H 

CM 

rH 

00 

rH 

oo 

CO 

rH 

tH 


in 

rH 

5 

r="J 

CD 

b 

O 

S.N. 

rH 

cq 

i 

• 

CO 

O 

rH 

CO 

co 

rH 

tH 

rH 

I 

LO 

0 

05 

rH 

1 

in 

CD 

« 

ci 

CO 

6 

IH 

LO 

1 

rH 

cq 

IH 

i 

cq 

CM 

1 


t- 

I 

© 

00 

c~ 


rH 

CD 

b 

IH 

rH 

rH 


tH 

in 

tH 

in 

in 

LO 

CO 

S' 

00 

r— 

05 

. _ | 

cq 


tH 


X 

cq 

05 

co 

tH 

co 

rH 

CD 

rH 

CO 

rH 

05 

CO 

CO 

LO 

LO 

o 

05 

i — 1 

00 

rad 

m 

£ 

£ 

CO 

• 

co 

1 — 1 
tH 

CD 

CD 

o 

cq 

rH 

• 

CD 

05 

rH 

CO 

00 

LO 

cq 

CD 

O 

rH 

• 

o 

cq 

cq 

t> 

co 

i 

CO 

rH 

tH 

cq 

CD 

00 

• 

05 

• 

rH 

CD 

« 

CD 


'<D 

b 

• 

cq 

1 

1 

i 

i 

1 


' 


v — " 

rH 

cq 

i 

| 



CM 



CO 

i 










1 


N " 




K 

§ 




















b 

o 




















u 

o 




















o 

CD 




















CO 

K 

CD 

b 

o 

05 


rH 


O 


LO 

S 

cq 


LO 

LO 

rH 

CD 

s 

CO 


cq 

cq 

cq 

CD 

tH 

tH 

(3.04 

CM 

rH 

CO 

co 


05 

cq 

in 

05 

rH 

05 


05 


Ph 

b 

o 

£ 

CD 

b 

o 

rH 

£ 

CD 

• 

cq 

cq 

05 

• 

CD 

1 

05 

05 

in 

t> 

rH 

CO 

in 

rH 

• 

in 

i 

LO 

• 

rH 

i 

CO 

• 

CD 

rH 

05 

• 

1 

00 

o 

CM 

• 

o 

LO 

IH 

CO 

• 

Si 

LO 

t 

rH 

rH 

05 

cq 

05 

• 

O 


CO 

i 








i 











b | 

a g 

f-l I u 

I 

<1 

b^ 


o 

rH 

o 


13 

o 

I— I 

CD 

b 

O 


m 


o 

rH 

t> 


b 

O 

£ 

1—1 

CD 


rH 


o 

o 

co 


a> cd 

b £t 
O & 

§ S 

£ 


co 

CD 

tH 






/-s 









. . 



tH 

rH 

co 

05 

OO 

rH 


rH 

00 

rH 

tH 

05 


CD 

cq 

05 

LO 

in 


© 

tH 

00 

CM 

CM 

in 

IH 

05 

rH 

05 

00 

in 

• 

CD 

cq 

© 

cq 

00 

© 

CO 

CM 

o 

o 

CM 

co 

© 

CO 

© 

00 

rH 

• 

© 

rH 

CM 

00 

rH 

rH 

• 

o 

l 

o 

\ 

rH 


CO 

rH 

rH 

1 

CD 

H 

co 

rH^ 

rH 

t 

rH 

si 

05 

1 

CM 

IH 

Si 



' — ' 



1 















oo 


cq 

o 

05 


o 

S 

CM 

IH 

© 

CM 

O 

rH 

LO 

CM 

CO 

o 

rH 

CO 

O 

CO 

rH 

LO 

cq 

00 

rH 

rH 

© 

rH 

© 

CM 

LO 

© 

o 

IH 

rH 

O 

CD 

cq 

CD 

00 

• 

© 

CM 

CO 

IH 

IH 

tH 

© 

rH 

00 

LO 

CM 

CD 

O 

O 

O 

CO 

cq 

CO 

£ 

© 

© 

CM 

© 

LO 

rH 

rH 

1 

rH 


o 

CM 

O 

© 

O 


o 


o 

co 



1 





1 








1 


CO 

0 



















rH 



o 


IH 

CM 

CO 

LO 

00 

r-“V 

00 

CM 

cq 

rH 

rH 


© 

© 

tH 

© 

% 

b^ 

o 

cq 

CM 

LO 

rH 

rH 

CO 

CO 

rH 

rH 

cq 

© 

cq 

CD 

rH 

tH 

IH 

00 

rH 

© 

£ 

rH 

00 

• 

CM 

• 

tH 

• 

rH 

CO 

• 

CM 

LO 

« 

05 

« 

© 

LO 

• 

LO 

rH 

CO 

© 

cq 

co 

• 

cq 

© 

r— < 

CD 

b 

o 

S.N 

o 

o 

CO 

rH 

tH 

Si 

05 

\ 

cq 

\ 

t 

rH 

\ 

rH 

1 

© 

CO 

rH 

rH 

1 

© 

1 

© 

© 


CM 


co 




LO 


co 


© 


146 




147 


Figure C5. — Lateral and Torsional Elastic Mode Shapes 
(Antisymmetric Boundary Conditions) 




MODE DEFLECTION (m/m) 


APPENDIX C - Continued 




(b) Lateral Component of Equivalent Booster 
Centerline Elastic Mode Shapes 


Figure C5. - Continued 
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APPENDIX C — Continued 




(c) Roll Component of Equivalent Booster 
Centerline Elastic Mode Shapes 


Figure C5. - Concluded 
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BOLL DEFLECTIONS (MODE 9) rad/m X 10' 



APPENDIX D 


SOLID MOTOR VEHICLE DATA 


This appendix contains a tabulation of the solid motor vehicle data that was used 

in this study. It also contains illustrations of the propellant sloshing tank geometries, 
the specific elastic vehicle models, and the elastic mode shapes. 

Essentially all of these data were generated during the Space Shuttle Phase E and 
follow-on studies. Therefore the data were developed using the customary (English) 
system of units. 

The definition of symbols is contained in Appendix A . 

The propellant sloshing parameters were generated by a version of the digital 
computer program described in Reference 2. 

The symmetric configuration uses a flexible beam model constrained to pitch- 
plane motion. The antisymmetric configuration also uses a flexible beam model, but 
this model has all of the antisymmetric degrees of freedom. The antisymmetric mode 
shapes were generated by the NASTRAN computer program. 



APPENDIX D — Continual 


TABLE Dl, - VEHICLE .MASS PROPERTIES 

Parameter Value (Liftoff + 10 seconds) 


T (3 orbiter engines) 

T S 

4. 893 X 10 6 N 
2.5 X 10 7 N 

(1. 1 x 10 6 lb) 
(5. 62 x 10 6 lb) 

M e (3 orbiter engines) 

10 200 kg 

(699 slugs) 

M r 

2.10 x 10 6 kg 

(1. 44 X 10® slugs) 

m t 

2.13 x 10 6 kg 

(1. 46 x 10® slugs) 

Ir (3 orbiter engines) 

25 260 kg-m 2 

(18 630 slug-ft 2 ) 

*xx 

6. 73 x 10 7 kg-m 2 

(4. 96 X 10 7 slug-ft 2 ) 

T yy 

3. 63 X 10® kg-m 2 

(2. 68 xlO® slug-ft 2 ) 

\z 

4. 12 x 10® kg-m 2 

(3. 04 X 10 8 slug-ft 2 ) 

^xz 

1.40 x 10 7 kg-m 2 

(1. 03 XlO? slug-ft 2 ) 


-23. 0m 

(-75.6 ft) 

^cS 

-23. 0m 

(-75.6 ft) 

^Ey 

6. 34 m 

(20. 8 ft) 

^Ez 

-7. 25 m 

(-23. 8 ft) 

% 

0. 94m 

(3. 08 ft) 

^RE 

1. 39m 

(4.56 ft) 

z cg 

-0. 49 m 

(-1.6 ft) 

5 op 

0. 09 rad 

(« 5 deg) 

5 oy 

0. 18 rad 

(« 10 deg) 

T RAT 

2/3 
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APPENDIX D - Continued 


TABLE D2. - STEADY-STATE TRAJECTORY DATA 


Parameter 

Value (Liftoff + 10 seconds) 

g 

9. 81 m/sec 2 

(32. 1? ft/sec 2 ) 

Q 

-0. 0089 rad/sec 

(-0. 51 deg/ sec) 

Uo 

44. 81 m/ sec 

(147 ft/sec) 

u 0 

14. 021 m/sec 2 

(46 ft/sec 2 ) 

v T 

44. 81 m/ sec 

(147 ft/sec) 

W o 

3. 92 m/ sec 

(12. 85 ft/sec) 

0 o 

1. 57 rad 

(90 deg) 

Mach Number 

0. 13 


Altitude 

2 72 . 8m 

(895 ft) 

Dynamic Pressure 

1245 N/m 2 

(26 lb/ft 2 ) 

TABLE D3. - 

■ TOTAL VEHICLE AERODYNAMIC DATA 

Parameter 

Value (Liftoff 

+ 10 seconds) 

x u 

-0. 000702 1/sec 



-0. 0010 l/sec 


Z 

u 

-0. 0024 l/sec 


z w 

0. 12 l/sec 


M u 

0. 000128 l/m-sec 

(0. 000039 l/ft-sec) 

M w 

0. 0021 l/m-sec 

(0. 00064 l/ft-sec) 

Y s 

-47. 55 m/sec 2 

(-156 ft/sec 2 ) 

L ,0 

-0. 0673 l/sec 2 


N 0 

-0. 028 l/sec 2 
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APPENDIX D - Continued 



1. ALL BULKHEADS ARE ELLIPTICAL SURFACES, 
MAJOR RADIUS = /2 X MINOR RADIUS 
2„ QUIESCENT FLUID LEVEL USED IN STUDY (~ -) 


Figure Dl. Main Propellant Tank Geometry 
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APPENDIX D — Continued 


TABLE D4„ - PROPELLANT SLOSHING DATA 


Propellant 

Tank 

Parameter 

Value 

Oxidizer 

M P 

18 607 kg 

(12 75 slugs) 

(L0 2 ) 

Lr, 

1. 01m 

(3. 3 ft) 


P 




^P 

16.6m 

(54. 5 ft) 


OJ p 

3. 75 rad/sec 

(0. 6 Hz) 


p 

1135, 7 kg/m 3 

(70.9 lb/ft 3 ) 

Fuel (LH 2 ) 

Mp 

5560 kg 

(381 slugs) 


L P 

2 . 06 m 

(6. 75 ft) 


A P 

1.93m 

(6. 33 ft) 


Wp 

2. 47 rad/sec 

(0. 39 Hz) 


P 

70. 16 kg/m 3 

(4. 38 lb/ft 3 ) 


TABLE D5. 

- AUTOPILOT AND ACTUATOR DATA 

Parameter 

Value 

K c 

12. 8 l/sec 

K o 

55 l/sec 

w c 

70 rad/sec 

w cn 

57 rad/sec 

£ cn 

0.48 

k r 

0. 7 sec 

k a 

variable 










TABLE D6.- ELASTIC VEHICLE MODAL DATA (PITCH PLANE BOUNDARY CONDITIONS) 
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Figure D4. - Elastic Model Schematic (Antisymmetric Boundary Conditions) 
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350.3 50.96 - 0.235 - 4.40 0.0886 0.297 9.09 0.000105 1.670 - 42.32 0.633 

( 24 ) ( 8 . 11 ) ( 0 . 0270 ) ( 0 . 0905 ) ( 0 . 509 ) (- 12 . 9 ) ( 0 . 193 ) 



TABLE D9. - YAW CONTROL PLANE AUTOPILOT SENSOR ELASTIC MODAL DATA 
(ANTISYMMETRIC BOUNDARY CONDITIONS) 


APPENDIX D - Continued 
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TABLE DIO. - ROLL CONTROL PLANE AUTOPILOT SENSOR ELASTIC MODAL DATA 
(ANTISYMMETRIC BOUNDARY CONDITIONS) 
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Figure D6. — Continued 







Mode 11, CO = 47. 69 rad/ sec (7. 59 Hz) 
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APPENDIX E 


THE ROOT LOCUS METHOD 


One of the most powerful tools applicable to the analysis and synthesis of launch 
vehicle flight control systems is the root locus technique. Its chief virtue is the graph- 
ical display of the motion of system roots as functions of various system parameters. 
This enables the control engineer to evaluate the effects of these changes directly . The 
root locations can be translated immediately into an estimate of the system time re- 
sponse or, if desired, its frequency response. The following discussion is an exposi- 
tion on the root locus method where the theoretical background is considerably abbrevi- 
ated. It can be found in varying states of clarity and completeness in almost any text- 
book on linear control theory; for example, References 4 and 5. 


Definitions 

The roots of the characteristic equation are a set of complex numbers that depend 
upon the coefficients of the equation. These roots are continuous functions of the co- 
efficients, hence of the system parameters. If the roots of the characteristic equation 
are plotted on the complex s-plane, they will trace out continuous curves as the several 
system parameters are varied. The path followed for a particular root is a root locus. 
For a given set of parameter values there is a single set of roots but through these 
roots a different set of loci will pass for each parameter variation. The question is 
then one of determining these loci. 


Theory 

The characteristic equation of a linear control system may be written in the form: 

m 

n (s - z.) 

i=l 

1 + GH (s) = 1 + K = 0 (El) 

n 

n (S - PH 
3=1 

where the z^ and pj are respectively the zeros and poles of the system, both real and 
complex, and K is the "high-frequency open-loop gain" — proportional to the low- 
frequency system gain. The function GH (s) is the open-loop transfer function for the 
system (see Figure El). Suppose the complex number, sq, is a root of Equation El. 
Then, since gain K is a real number, the characteristic equation is equivalent to two 
relations: 
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C(t) 

Controlled Variable 

R(t) 

Reference (Command) Variable 

E(t) 

Error Variable 

B (t) 

Feedback Variable 

G(s) 

Forward Loop Transfer Function 
(Including Compensation) 

H (s) 

Feedback Loop Transfer Function 
(Including Compensation) 


Figure El. Generalized Block Diagram for Control Systems 


K 



(E2) 


Arg 


m 

n (s 0 - z ) 
1=1 


n 




^Arg^-z.) "S Arg(s 0 -P j )=(2r + I)j 


(E3) 


j=l 
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APPENDIX E - Continued 


where r is any integer. If a point can be found in the s-plane satisfying Equation E3, 
there will be only one value of K for which Equation E2 is satisfied. The root locus for 
Sq under variation of the open-loop gain is the locus of points satisfying Equation E3 

with K a parameter. 

This procedure will define the motion of the system roots as a function of K. A 
few observations are in order. On the s-plane there may exist "saddle points, " points 
where the derivative of the complex gain function given by 

n 

n (s - p ) 

tl 

m 

n (s - z.) 

i=l 




goes to zero. If the locus passes through a region surrounding a saddle point it will be 
found that the distance along the locus for a given percentage change in gain increases 
dramatically; furthermore, in plotting the locus the points s Q lying on the locus will be 
difficult to find. This situation usually occurs when two separate loci approach each 
other from opposite directions, when it will be difficult to determine if the locus is as 
shown (for example) in Figure E2a or E2b. From the engineering point of view, it makes 
little difference; the location of a root lying in the region of a saddle point will be 
sensitive to parameter variations, and the system should be designed accordingly, 
i. e. , to be insensitive to the precise location of the root. This figure shows a typical 
situation encountered when analyzing the low-frequency parasitic modes (represented 
by the pole-zero dipole near the jeo axis). If the operating gain is either K-^ or K^, the 
situation is clear. But an intermediate value of gain will result in two complex con- 
jugate pairs of roots within the region of the saddle point that are difficult to locate 
with precision. This can best be appreciated by one familiar with the use of the 
Spirule: the locus is hard to "pin down. " The difficulty is less severe if the brute- 
foi’ce methods of machine computation are used, but even here the number of iterations 
required to locate the roots will be somewhat larger than usual. It is not at all unusual 
to run across difficult cases where sophisticated computational techniques are required; 
the ordinary routines fail to converge. 

For synthesis purposes, one is concerned with the effect of lead or lag compen- 
sation in the open-loop transfer function on the root location or, what amounts to the 
same thing, on the configuration of the loci. In essence the effect is as follows: in the 
upper half of the s-plane, loci emanating from complex poles will rotate clockwise 
about the pole as lag compensation is added. They also will tend to decrease in length 
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a. Variation 1 b. Variation 2 


Figure E2. Variations in Root Loci Behavior in the Region of a Saddle Point 

(the closed-loop roots move closer to the open-loop poles) by an amount proportional 
to the attenuation introduced by the lag. The amount of rotation is (approximately) 
equal to the phase lag introduced at the frequency corresponding to the system root of 
interest. In terms of Figure E2, the situation shown in (a) represents somewhat more 
lag than (b). It also shows the effects of lag on loci emanating from poles on the real 
axis: the loci will tend to move toward the right-half plane. For lead compensation, 
the situation is exactly reversed: the loci tend to increase in length and rotate counter- 
clockwise. 

Phase and gain margins are conveniently measured for parasitic mode roots on 
root locus plots if the definitions in Figure E3 are used. One draws a vector between 
the pole and the nominal closed-loop root location, and inscribes a circular arc about 
the open-loop pole of radius equivalent to the distance of the closed-loop root from the 
pole. A vector is drawn from the pole to the intersection of the arc with the jto axis. 

The angle 0 m between the two vectors is the phase margin. The gain margin is the 
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Figure E3. Gain and Phase Margin, Root Locus Interpretation 

ratio between the gain at the point where the locus crosses the axis and the gain at the 
operating point (in decibels): 


*2 

Gain Margin = 20 log 1Q (E5) 


Time and Frequency Response 

The root locus technique is fundamentally a means of locating the roots of the 
characteristic equation. Its chief virtue is the display of these roots as functions of 
various system parameters. But most important are the system response properties 
to command or disturbance inputs, both inferred from the location of the open-loop 
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APPENDIX E - Continued 

system poles and. zeros. The s-plane diagram of these poles and zeros serves as a 
graphical means for determining the time and frequency response of the system. 

The transfer function of a system can be represented by 

m 

K n ( S - z.) 

T ( S ) = — (E6) 

n 

n (s - p i ) 

3=1 

The input can be a command, a disturbance, or a test input; the output is the variable 
whose response is of interest. Consider the frequency response obtained by substi- 
tuting jco for s, 


T (jco) 


m 

K n (jco - z.) 
i=l 
n 

n (jw - p.) 

1=1 


T (jco) | 


j ArgT(jto) 

c 


(E7) 


This may be evaluated graphically from the s-plane plot of the poles and zeros for any 
desired to by picking a point on the positive jco axis corresponding to the desired fre- 
quency and measuring off distances to all the poles and zeros and their respective 
phase angles according to 


|T (jco), n 

n |jw - pJ 
3=1 J 

m n 

ArgT (jco) = 23 Arg “ z i) " 23 Arg " p j> (E8) 

i=l j=l 


m 


K n |jw - z.| 
i=l 
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where the angles are measured at the poles and zeros counterclockwise from the 
3 o’clock position to the j to point on the jto axis. 

Consider the transform of the time response to an input, C(s). 


C (s) = R(s)T(s) 


m 

K n <s - z.) 

— R (s) (E9) 

n (s - pj) 
j=i 


If R (s) is in closed form having only first-order poles (this excludes ramp inputs, for 
example) at frequencies different from the system poles, the response may be evalu- 
ated from an s-plane plot of the poles and zeros of C (s). C (s) may then be represented 
by: 


k 

k' n (s - z.) 

i=l 

C (s) = 1 

n (s - p i ) 

j=i 


(E10) 


where the indices k and A (rather than m and n) indicate additional poles and zeros in 
C (s) due to R (s). C (s) has the partial fraction expansion, 


C(s) 




(s - Pj) 


(Ell) 


where the 



residue is 



(s - p p ) C (s) 


s=p 


P 



(E12) 
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That is, 

|a| = K'n|p - z.|/ n |p - Pj| 
i=1 j=i 

_ j/p 

k i 

0 = E Ar § (p P - z i> - £ Ar & (p p - Pj) 
i=l j=l 

j/p 

In short, the residue at the p™ pole of C (s) is determined by measuring the distance 
and angles as before from the s-plane plot. The time response is given by 

Pjt 

C (t) = V A.e 3 (E14) 

j=l 

If some of the poles are complex, the residues appear in complex conjugate pairs. 
Suppose there are poles at p 1 = + j^) and p 2 = (0^ - jo^) with residues A ± and 

Ai e - 3^. The corresponding contribution to the time response C (t) is 

C-^ (t) = 2 A^ e 1 cos (co^ t + 0) (E15) 

The corresponding time functions may be graphically added to give the overall time 
response. 
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